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LINEARLY CONSTRAINED EVOLUTIONS OF CRITICAL POINTS AND 
AN APPLICATION TO COHESIVE FRACTURES 

M. ARTINA, F. CAGNETTI, M. FORNASIER, F. SOLOMBRINO 


Abstract. We introduce a novel constructive approach to define time evolution of critical 
points of an energy functional. Our procedure, which is different from other more established 
approaches based on viscosity approximations in infinite dimension, is prone to efficient and con¬ 
sistent numerical implementations, and allows for an existence proof under very general assump¬ 
tions. We consider in particular rather nonsmooth and nonconvex energy functionals, provided 
the domain of the energy is finite dimensional. Nevertheless, in the infinite dimensional case 
study of a cohesive fracture model, we prove a consistency theorem of a discrete-to-continuum 
limit. We show that a quasistatic evolution can be indeed recovered as a limit of evolutions of 
critical points of finite dimensional discretizations of the energy, constructed according to our 
scheme. To illustrate the results, we provide several numerical experiments both in one and two 
dimensions. These agree with the crack initiation criterion, which states that a fracture appears 
only when the stress overcomes a certain threshold, depending on the material. 
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1. Introduction 

In this paper we introduce a novel model of time evolution of physical systems through linearly 
constrained critical points of the energy functional. In order to include all the applications we 
have in mind, we consider both dissipative and nondissipative systems. 

Our approach is constructive, and it can be numerically implemented, as we show with an 
application to cohesive fracture evolution. Since we are eventually interested in being able to 
perform reliable numerical simulations, we consider at first finite dimensional systems. We then 
also give a concrete case study showing how our results can be adapted to describe infinite 
dimensional systems. 

Below we recall the general framework, to which we intend to contribute. Then we present and 
comment our results, also in comparison with related contributions appeared in recent literature. 

1.1. Critical points evolutions in the literature. When describing the behaviour of a physi¬ 
cal system, one can try may want to describe it through the time evolution of absolute minimizers 
of the energy. This modeling has been pursued, for instance, in laiiiiiiaiisiiiTiiiEiiiH]. From 
an abstract point of view, this amounts to requiring that a global stability condition is satisfied 
at every time. In this case, the notion of solution fits into the general scheme of energetic so¬ 
lutions to rate-independent systems (see [26]). However, it is not always realistic to expect the 
energy to be actually minimized at every fixed time. In fact, global minimization may lead the 
system to change instantaneously in a very drastic way, and this is something which is not very 
often observed in nature. For this reason, several authors recently introduced time evolutions 
only satisfying a local stability condition for the energy functional HlElElliaiMlEIllMlESI. 
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More precisely, given a time-dependent functional F :Y x [0,r] —)• M, where y is a Banach 
space, an evolution of critical points of F is a function u : [0,T] ^ Y which satisfies 

0 G duF{u{t),t), for a.e. t G [0,T], (1-1) 

where duF denotes the subdifferential of F with respect to u. Typically, the existence of such 
an evolution is proven by a singular perturbation method. More precisely one first considers, 
for every e > 0, the e-gradient flow 

— G duF{u^{t),t) (1.2) 

with an initial datum u^(0) = uq, where uq is a critical point of F(-,0). Then, passing to the 
limit as e —>• O"*", converges to a function u satisfying 

We now give a detailed description of our results and we comment on them. 


1.2. Setting of the problem and main results. In this paper we consider the evolution of 
a system which is driven by a linear external constraint. This can model several of situations 
of interest, such as prescribed boundary data, integral constraints, or the coupling with a linear 
(partial) differential equation. We will state and discuss the problem in a discrete (finite¬ 
dimensional) setting, where our main results are obtained. Later on we will also comment 
on how to possibly recover solutions to problems defined in infinite dimension. 

Let £ and F be two Euclidean spaces of dimension n and m, respectively, with n > m. 
We want to study the evolution in a time interval [0, T] of a physical system, whose states are 
described by vectors u G T, and whose energy is given by a function J : £ —)■ [0,-|-oo). We 
assume that a time dependent constraint / : [0, T] —?■ is imposed. More precisely, the only 
admissible states at each time t G [0, T] satisfy Av = f{t), where A : £ ^ F is a surjective linear 
operator. In addition, we consider a convex positively 1-homogeneous dissipation functional 
ij) '. £ ^ [0, oo). In order to consider at the same time dissipative and nondissipative systems, 
we introduce a switching parameter a G {0,1}. 

Before stating the main results of the paper, let us give introduce some notions, which are 
useful in the rest of the paper. At any time t G [0, T] and for a fixed v £ £, one can consider 
the problem 

min {J{v) + a^p{v — v)}. (1-3) 

Av=f(t) 

If u G T is a minimizer for ( |I.3[ ), then 

Au = f{t) and {dJ{u) + d{aip){u— v)) r\vaTi{A*) ^ 0, (1.4) 


where A* ■. F ^ £ is the adjoint of A, ran(A*) denotes the range of A*, and 
subdifferential of J at u. A critical point of v ^ J{v) + a'ip{v — v) on the 
is any vector u £ £ satisfying (1.4) where, for every f £ F, we set 


dJ{u) is the Frechet 
affine space A{f{t)) 


A(/) := {v £ £ : Av = /}. 


A discrete quasistatic evolution with time step 6 £ (0,1), initial condition vq £ T, and constraint 
/, is a right-continuous function vs '■ [—6,T] -£ £ such that 

• vsit) = Uq for every t £ [—5, h); 

• vs is constant in [0, T] H [f<5, [i + l)h) for all z G Nq with i6 < T; 

• vs{i6) is a critical point of u i—;■ J -|- afj{v — vs{{i — 1)<5)) on the affine space A{f{i6)) for 
every i G No with i5 < T. 
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Moreover, we say that a measurable function u : [0, T] —)• £1 is an approximable quasistatic 
evolution with initial condition vq and constraint /, if for every t G [0, T] there exists a sequence 
0 ^ and a sequence {vsi_)k£n of discrete quasistatic evolutions with time step 5^, initial 
condition vq, and constraint /, such that 


lim inf \vs^^{t) - v{t)\g = 0 . 

fc^+OO 


We are now ready to state the main results of the paper (see Theorem 2.15 and Theorem 2.16). 

Dissipative systems. Let vq be a critical point of u i—)• J{v) + 'ip{v — vq) on the affine space 
A(/(0)). Under suitable assumptions on J, if;, A, and / (see Theorem 2.15| ) we prove that there 
exist V G BV{[0,T]-,£) of bounded variation and q G L°°([0,T]-,iF') such that: 

(A) v{-) is an approximable quasistatic evolution with initial condition vq and constraint /; 

(B) Local stability: A* q{t) G dJ{v{t)) + 9V’(0) for C^-a.e. t G [0, T] 

(C) Energy inequality: the function s i—>• {q{s), f (s)) jr belongs to L^(0,T) and 

J{v{t2)) +yaT^{v;[ti,t2]) < J{v{ti)) + f {q{s),f{s))jrds, 


for every 0 < ti < t 2 < T, where the ?/;-variation Var. 0 (u; [ti,t 2 ]) is defined by (2.9) and 
(•, •)j- denotes the duality product in F. 

Any function v G BV{\t),T];£) satisfying (A), (B), and (C) is a rate independent evolution. 
We also remark that, as a consequence of the local stability (B), the energy inequality becomes 
an equality in all the nontrivial intervals where the solution v{t) is absolutely continuous (see 
Theorem 2.18, where also the nondissipative case is treated). In addition, in such intervals the 
doubly nonlinear inclusion A*q{t) G dJ{v{t)) + ^^l^{v{t)) is satisfied. It is however a well-known 
fact that, due to nonconvexity, our solutions can in general develop time discontinuities, where 
additional dissipation appears (see m)- 

Non dissipative systems. When there is no dissipation we can still prove an existence result, 
although the evolution obtained is in general expected to be less regular in time. This is due to 
the fact that the absence of dissipation causes loss of compactness, since simple estimates of the 
total variation of the approximating solutions are now missing. This is undoubtedly a point of 
great interest in the analysis of the degenerate case (see also [ 2 ] for an abstract approach in the 
unconstrained setting). To compensate the loss of compactness, we need to make an additional 
assumption on the energy functional J (see condition (J4) in Section]^. If vq is a critical point 
of J on the affine space A(/(0)), we prove (see Theorem 2.16) that there exists bounded and 
measurable functions v : [0, T] —)■ T and (7 : [0, T] —)• T"' such that: 

(a) v{-) is an approximable quasistatic evolution with initial condition vq and constraint /; 

(b) Local stability: A*q{t) G dJ{v{t)) for every t G [0,T]; 

(c) Energy inequality: the function s 1 —>■ {q{s), f (s)) jr belongs to L^(0,T) and 

/*i 


J{v{t))<J{vo)+ {q{s),f{s))jrds, 


(1.5) 


for every t G [0, T]. 


^More precisely, by constructing suitable representatives of q{t), one can even ensure that local stability holds 
at all continuity points of t 1 —>■ v{t). We do however not deal in detail with this technical aspect. 
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We call any measurable and bounded function u : [0,T] —>■ satisfying (a), (b), and (c) a 
weak potential type evolution. Evolutions of this kind (although without constraints) have been 
widely considered in literature, as limits for e —)• 0 of gradient flows of the type as in |39j ' 
or of systems with vanishing inertia (see mm)- We observe that the term {q{t), f (t)) jr in (c) 
physically corresponds to the virtual power due to the external constraint. In the case where J 
were smooth, thanks to (b) this term could indeed be rewritten as {DJ{u{t)),w{t))j^, for any 
smooth curve w{t) with Aw{t) = f{t). 

Note also that in our definition the precise value of v at every point t matters. In particular, 
the initial condition u(0) = vq has a meaning, and the energy inequalities (C) and (c) need to be 
satisfied at every time. Thus, we are in general not identifying functions differing on null sets, 
as it is usual in U’ spaces. 

We point out that the main novelty in our approach is not in the existence results per se, 
but rather in the constructive algorithmic procedure that we provide, see (3.3). Notice that, 
differently from the vanishing viscosity approach, the parameter rj appearing in (3.3) remains 
fixed, throughout all the algorithm. Heuristically, our inner loop aims at finding the nearest 
critical point through a sort of discretized instantaneous generalized gradient flow, see Figure 
and Figure This can be obtained by looking at the long time behavior of the minimizing 
movements of the functional, for a fixed time step 5. More details about this point are given in 
the next subsection. 



Figure 1. A pictorial idea of the algorithm when a = 0. Let z G N be fixed, 
and suppose that is a critical point of J on the affine space A(/((z — 1)(5)). 
The algorithm allows to pass from to a vector u*, which is a critical point 

of J on the affine space A(/(z5)). With abuse of graphical notation, constrained 
critical points are identified in the picture with those points where the slope of 
the functional J is zero. A magnified version of the area delimited by the blue 
rectangle can be seen in Figure 


As explained below, we also show how it is possible to apply our results to a model of cohesive 
fracture evolution. In this specific application, we decided to work in the nondissipative setting, 
for several reasons. First of all, many results are available in the literature in this case (see 
0), and this allows us to test our methods. Secondly, the notions of fracture evolution which 
are used to model the cohesive dissipative case are not easy to use in applications, since they 
require rather delicate tools of functional analysis (a formulation is spaces of Young measures, 
see [7]). Finally, this is a relevant application for the degenerate case, where BV estimates are 
not available. 
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Figure 2. A pictorial description of the sequence provided by the al¬ 

gorithm in (3.3). For each j G N, the vector Vj is the unique minimizer of the 
strongly convex function Jn^vj^iiy) '■= J{v) + r/|u — When j —)> oo, up to 

subsequences, we have that Vj —)■ u*, where u* is a critical point of J in the affine 
space A{f{i6)). 


We are aware that a more realistic model should take account of the monotonicity of the crack 
growth, as done for instance in [3 da EH]. This issue can be dealt with in the case of brittle 
fracture, thanks to the Jump Transfer Lemma. Unfortunately, to date, this tool is not available 
in the cohesive case (see |18] for details), and even if it were existing, it would not help with the 
lack oi BV estimates. 


1.3. Comments on our result. We would like to emphasize a few relevant and defining aspects 
of our results. 


(i) We prove the existence of an approximable quasistatic evolution for a large class of energy 
functions J and linear operators A (see (2.1) and conditions (J1)-(J4) in Section]^. In 
particular, these include nonsmooth and nonconvex energy functionals. 


(ii) We stress the fact that v{t) is supposed to visit at different times t critical points of the 
functional v i—>■ J{v) + a^l^{v — v{t)) over the affine space A(/(t)). This condition is rather 
general, if compared to the usual requirement of focusing on global minimizers of J over 
A(/(t)). An evolution along critical points is in general more realistic and physically 
sound. Moreover, it is important to notice that the viscous approximation usually does 
not provide an evolution along critical points, unless one lets the viscosity to vanish. 
Hence, in contrast to this well-established approximation method, our approach starts 
with a truly consistent approximation from the very beginning. However, the analysis 
of such an evolution is usually very involved and, in absence of dissipation, may allow 
for solutions v that are just measurable in time. While we shall be content with the 
generality of our results, we have to live with the fact that our solutions may not be 
regular in the nondissipative case. 


(iii) Our approach is constructive. This is an important fact since, as we shall emphasize 
later, the functional v i—)■ J{v) + aip{v — v{t)) may have multiple feasible critical points 
at the same time t. Hence, in order to promote uniqueness of evolutions, or even just 
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their measurability, we need to design a proper selection principle. Accordingly, we 
shall design our selection algorithm in such a way that the selected critical point is the 
closest - in terms of Euclidean distance - to the one chosen at the previous instant of 
time, unless it is energetically convenient to perform a “jump” to another significantly 
different phase of the system (see formula (3.3)). This corresponds to a rather common 
and well-established behavior of several physical (and non physical) systems (see, for 
instance, [6l Section 9]). Moreover, this method can be easily implemented by means 
of a corresponding numerical method [3|, as we will show with the above mentioned 
application to cohesive fracture evolution, see Section Thus, our evolution is the 
result of a constructive machinery (an algorithm), which is designed to emulate physical 
principles, according to which a critical point is selected in terms of a balance between 
neighborliness (accounting the Euclidean path length between critical points) and energy 
convenience. In our view, this feature is of great relevance as we provide a black box, 
whose outputs are solutions. 

(iv) As already mentioned, the proof of the main results (Theorems 2.15 and 2.16) is given in 
a finite dimensional setting. This is due to the fact that, in the inhnite dimensional case, 
the subdifferential is in general not closed with respect to the weak convergence in the 
domain of the energy. Such a difficulty could be overcome by requiring that the energy 
functional has compact sublevels, an assumption which is quite common in literature, 
provided the domain of the energy is suitably chosen. On the other hand, the choice 
of a weaker topology for the domain may not always comply well with other conditions 
on the energy, that we need to prove the existence results. This is in particular the 
case of the key condition (J3) (see Section]^, which allows to control the virtual power, 
due to the external constraint, in terms of the energy. While not representing a major 
hurdle in the uncostrained case with time-dependent energy functionals, this issue seems 
particularly relevant for the problem we study (see Remark 6.10 for further details). This 
motivates our choice of first dealing with a finite dimensional setting, and then extending 
the results to our model case (see Theorem 7.5) with a problem-specific technique. 


(v) As an important remark, we stress that in general all of the constants appearing in 
the technical assumptions in Section could depend on the dimension of the consid¬ 
ered Euclidean spaces. Thus, our results can be applied to physical systems that can 
assume (a discrete or a continuum of) infinitely many states, provided all the relevant 
estimates obtained are dimension free. Eor this reason, we state very clearly which are 
the parameters affecting the constants that come into play in the crucial proofs (see 
Remark |4.3| and Remark |5.2[ ) . We give an important application to cohesive fracture 
evolution (see Section and Section , showing how also infinite dimensional systems 
can be approached with our method. In particular, we eventually provide an alternative 
proof of the existence of evolutions of critical points for the cohesive fracture model 
firstly proven in [6]. 


(vi) The numerical simulations that we provide in Section for the cohesive fracture evolu¬ 
tions agree with physically relevant requirements, such as the crack initiation criterion 
(see 0 Theorem 4.6]), which states that a crack appears only when the maximum sus¬ 
tainable stress of the material is reached. We also mention that numerical simulations, 
obtained instead with the vanishing viscosity approach, have already appeared in [24] . 
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To the best of our knowledge, this is the first time that an algorithm providing critical points 
evolution is introduced in such a generality, especially to treat consistently and stably also 
nondissipative models. 

However, it is worth mentioning that our approach, although derived independently, resembles 
similar methods which have appeared recently in the literature. In [32], a related scheme has 
been for instance investigated in order to obtain a general existence result in a nonconvex but 
smooth setting. The author also takes into account viscous dissipation effects, and provides a 
constructive time rescaling, where the evolutions have a continuous dependence on time. This 
idea, in particular, generalises previous approaches for systems driven by nonconvex energy 
functionals (linillllETlETI). Moreover, the author shows an approximation result that in spirit 
is close to our Theorem 7.5, and to previous results in |29l[36|. However, the results in [32| are 
obtained under the assumption of regularity of the energy functional, and in an unconstrained 
setting. In particular, stability of critical points after passing to the limit is recovered through 
a very strong assumption f |32l ( 8 ), Theorem 2.3]), which would seem quite unnatural for a 
constrained nonsmooth problem. 

Concerning other contributions, we also mention that a very general incremental minimization 
scheme, involving a quadratic correction with a fixed parameter ^ > 0 , has been just proposed 
in |30j . even in connection with abstract dissipation distances in metric spaces. 


Another algorithm, showing some analogies to (3.3), has furthermore been recently considered 
for a case study of phase field fracture coupled with damage in |22j . In this case, the energy 
is nonconvex, but separately convex in the two state variables. Therefore, instead of adding a 
regularization, the authors define the discretized evolution through fixed points of an alternate 
minimization scheme. Also in this case a time reparametrization, where a full energy-dissipation 
balance holds, is provided. The exploitation of similar techniques in connection with our problem 
is another interesting issue that we plan to pursue in the future. 

The plan of the paper is the following. In Section]^ we state the main results of the paper. 
Theorem 2.15 and Theorem 2.16, whose proofs are given in Sectionj^and Section]^ respectively. 
Section [6 is devoted to the description of the cohesive fracture model introduced in |B] . In the 
same section we introduce a space mesh and spatially discretize the problem. In Section]^ we 
pass to the limit as the size of the mesh tends to 0 , thus obtaining a new proof of the result in 
| 6 l Theorem 4.4]. Finally, numerical simulations are given in Section]^ 


2. Setting of the problem and main result 

2.1. Basic notation. Throughout all the paper, we use the notation No ;= N U {0}, and we 
denote by the standard Lebesgue 1-dimensional measure in M . Let £ ~ M"' and T ~ 
be two Euclidean spaces with dimension n and m, respectively, with m < n. We consider an 
energy function J : T —>■ [0, -boo), a linear operator A : T —)• and a time dependent constraint 

/ : [0, T] —>■ iF. We will assume that A is surjective. Equivalently, we will suppose that that 
there exists 7 > 0 such that the adjoint operator A* : F' ^ £' satisfies 

\A*(1\£' > iIqIt' for every q € F'. (2.1) 

We will frequently use the space BV{[a,b], X) of functions of bounded variation from a time 
interval [a, b] to a Banach space X. 

Before proceeding, let us recall some basic notions of differential calculus that are used in the 
sequel. Given u G £ and S : £ ^ [0,oo), we recall that the Frechet subdifferential dS{u) C £' 
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of 5 at a point u £ S is defined in the following way: 

S{v) - {S{u) + {C,v- u)s) 


e G dS{u) 


0 < lim inf 

V^U 


\v - u\e 

where {-,■)£ denotes the standard scalar product in S. For every f £ T, we set 

A(/) ■.= {v££:Av = f}. 

We can now give a precise definition of critical point in our setting. 

Definition 2.1. Let f £ T and S' : <5 —>■ [0,oo). Assume that S has nonempty subdifferential 
at every point. We say that u £ £ is a, critical point of S on the affine space A(/) if 


Au = f and dS{u) n ran(A*) ^ 0, 
where ran(A*) denotes the range of A*. 


( 2 . 2 ) 


Remark 2.2. One can check that condition (2.2) implies, in turn, that 

^ ^ T • r r 1 / 4\ 

U < liminf ■ ■ —— ^ tor every w £ ker(A), 

£—>• 0 + £ 

where ker(A) denotes the kernel of A. 

Remark 2.3. The definition above is motivated by the fact that if u £ £ satisfies 

min S{v) = S{u), 

Av=f 


(2.3) 


then (2.2) holds true. 

We also recall that a function S : £ —)• [0, oo) is said to be strongly convex if there exists v > 0 
such that whenever , t :2 £ £ we have 

u\vi - n 2 || < (Cl - - V 2 )£ for every Ci G dS{vi) and C 2 G dS{v 2 ). 


2.2. The energy fnnctional. We can now state our assumptions on the energy functional J. 
We suppose that: 

(Jl) the functional v 1 —> J{v) + |Ai;|3r is coercive; 

(J2) there exists t] > 0 such that v 1 — > Jr]{v) := J{v) + r]\v\g is strongly convex; 

(J3) there exists L > 0 such that, for every v £ £, £ dJ{v) => ICI^:' < L (J(n) + 1). 

Before proceeding with the setting of the problem, we make some remarks on the assumptions 
above. 

Remark 2.4. Condition (J2) implies that J is a smooth perturbation of a convex function and, 
therefore, it is locally Lipschitz continuous. From this, it follows that J is almost everywhere 
differentiable and its Frechet subdifferential is nonempty at every point. 

Remark 2.5. Let 7 / > 0 be such that (J2) holds true, and let v £ £. Then, the functional 
Jn^v : T —[0, 00 ) defined as 

Jrj,v{v) ■= J{v) + r]\v — Till for every v £ £, (2-4) 

is also strongly convex, with the same constant v. 
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Remark 2.6. If J is globally Lipschitz continuous with Lipschitz constant L, then 

^ G dJ{v) < L, (2.5) 

for every u G <5. Moreover, if J = Ji + J 2 where Ji is lower semicontinuous and J 2 is of class 
then the decomposition 

dJ{v) = dJi{v) + DJ 2 {v) (2.6) 

holds true at every point v £ £ such that dJ{v) / 0, where DJ 2 {v) denotes the Frechet 
derivative of J 2 at u. 

Remark 2.7. As shown in [3l Remark 2.5], it suffices to check condition (J3) only at differen¬ 
tiability points of J, to ensure that it is satisfied at every point. 

2.3. Dissipation Functional. We introduce now the dissipation functional if -. £ ^ [ 0 , 00 ), 
which measures the energy, which is lost when passing from one state to another one. We assume 
that if satisfies the following: 

('hi) if{v) = 0 if and only if u = 0 ; 

('h2) il){vi + V 2 ) < if{vi) + if{v 2 ); 

('I'3) ip{Xv) = Xif{v) for every A > 0 and v £ £; 


Remark 2.8. Note that assumptions ('h2) and ('I'3) imply that if is convex. Note also that in 
general if is not a norm, unless symmetry holds (i.e., if{—v) = if{v) for every v £ £). 


From the assumptions above it follows that there exists a constant c > 0 such that 


1 


-j'l’k < ' 0 (^) ^ cjujf for every v £ £. 


(2.7) 


and the references therein) 

( 2 . 8 ) 


Under these assumptions, it is well known (see, for instance, 
that 

0 / dif{v) C K* for every v £ £, 

where we set 

K* := dif{0) = {(,££': {f,,w)£ < if{w) for every w £ £}. 

One can also check that, for every v £ £, the subdifferential dif{v) of if at v is characterised by 

dif{v) = K* n{f,£ £' : (^, v)£ = if{v)}. 

Remark 2.9. The set K* is convex (as it is a subdifferential) and bounded; indeed, the inclusion 
K* C R(0,c) (this one being the ball with radius c in £') simply follows from (Tl) and (2.7). 


Remark 2.10. Since 'h is proper convex, from (J2), the convexity of (2.6), and the Moreau- 
Rockafellar Theorem f |35l Theorem 23.8]), the decomposition 

d{J -|- if){v) = dJ{v) + dif{v) 

holds at every v £ £. 


For if as above and u £ BV{[a, h],£), the ?/:-variation of u is defined as 

Var^(u; [a,b]) = sup (V’(ri(U)) “ 'if{u{ti-i))) : a = to < U < • • • < 4 = &, fc G n| . (2.9) 

If one takes i/: = j • j in the above definition, one retrieves the usual definition for the pointwise 
variation of a function. For all a < c < 6 the equality 

Var^(u; [a,b]) = Var^(u; [a,c]) Var. 0 (u; [c, 6 ]) 


( 2 . 10 ) 
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immediately follows from the definition and the subadditivity oiil). If rt is additionally absolutely 
continuous, it is well known that 

Var^(ti; [a, 6]) = f '4>{u{s))ds. (2-11) 

J a 


2.4. Main results. In the following, together with rate independent evolutions, in which dis¬ 
sipation is present, we will also consider weak potential type evolutions, where there is no dis¬ 
sipation. From the technical point of view, the absence of dissipation translates into a lack of 
compactness. For this reason, we need an additional assumption to treat this case. We shall 
assume 


(J4) There exists a positive constant Cj^r^ > 0 such that for every u G T 
A*qi G dJrj^vivi), i = l,2 {qi - q2, Avi - Av2)j^ < Cj^rj\vi - V2\e\Avi - Av2\j^. (2.12) 


Remark 2.11. Although condition (J4) above might seem quite technical, it is automatically 
satisfied when J G C^’^. Indeed, in this case dJ{v) is single valued at every v G £, and coincides 
with the differential DJ{v). Then, denoting by M the Lipschitz constant of DJ{-) and using 
(2.1), one has 

\qi - g2|j-/ < -\A*qi - A*q 2 \s' = -\DJ{vi) - DJ{v 2 )\s' < —\vi - V 2 \s- 
7 7 7 


At this point, (2.12) simply follows by the Cauchy-Schwarz inequality. 

Remark 2.12. We will sh ow in a concrete example that condition (J4) can also be satisfied 
when J ^ (see Section 6.3). 


Before stating our main results, we give again and in more detail the notion of discrete and 
approximable quasistatic evolution, respectively. When this is possible, in the following we treat 
at the same time the cases with and without dissipation. To this aim, we introduce a switching 
parameter a G {0,1}, in such a way that a = 0 corresponds to the situation without dissipation, 
while in the case a = 1 dissipation is present. 


Definition 2.13. Let a G {0,1}, let uq G T be a critical point of J -|- — uq) on the affine 

space A(/(0)), and let 5 > 0. A discrete quasistatic evolution with time step 6, initial condition 
Vo, and constraint / is a right-continuous function vs '■ [—6,T] ^ £ such that 

• vsit) = vq for every t G [—5, 5); 

• Vs is constant in [0, T] (7 [i<5, [i + 1)<5) for all i G Nq with i6 < T; 

• vs{i6) is a critical point of u i—)■ J -|- aip{v — vs{{i — 1)<5)) on the affine space A{f{i6)) for 
every i G Nq with i5 < T. 


Definition 2.14. Let a G {0,1} and let uq G T be a critical point of u i—)• J(v) -G atplv — vq) 
on the affine space A(/(0)). A bounded measurable function v : [0,T] ^ £ is said to be an 
approximable quasistatic evolution with initial condition vq and constraint /, if there exists a 
sequence —)• 0^ and a sequence {vs^{t))k^^ of discrete quasistatic evolutions with time step 
5k, initial condition vq, and constraint /, such that, for every t G [0, T], 

lim inf \vs^it) - v{t)\£ = 0. (2.13) 

fc^+OO 

We are now ready to state our main results. The first one is an existence result for rate 
independent evolutions. 
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Theorem 2.15 (Existence of rate independent evolutions). Let a = 1, and suppose that 
(2.1), (Jl), (J2), and (J3) are fulfilled, and that (^1), ('if2), and (if3) hold true. Let f G 
W^{[Q,T];F), and let vq be a critical point o/u i—)• J + fi{v — uq) in the affine space A(/(0)). 
Then, there exist v G BV{[Q,T]]£) and q G L°°{[0,T]-,F') such that: 

(A) v{-) is an approximable quasistatic evolution with initial condition vq and constraint f; 

(B) A*q{t) G dJ{v{t)) + K* for C^-a.e. tG[0,T]; 

(C) The function s i—)■ {q{s), f (s)) jr belongs to L^{0,T), and for every 0 <ti <t 2 <T 

J{v{t 2 )) +yax^{v;[ti,t 2 ]) < J{v{ti)) + f {q{s), f{s))jr ds. 

Jti 

In the case without dissipation we need to add the additional assumption (J4), and we obtain 
measurability, but in general no further regularity, of the evolution. 

Theorem 2.16 (Existence of weak potential type evolutions). Let a = 0, and suppose that 
(2.1), (Jl), (J2), (J3), and (Jf) are satisfied. Let f G W^’‘^{[d,T\,F), and let vq be a critical 
point of J in the affine space A(/(0)). Then, there exist bounded and measurable functions 
u ; [0, T] —)■ and q : [0, T] —)■ F' such that: 

(a) v{-) is an approximable quasistatic evolution with initial condition vq and constraint f; 

(b) A*q{t) G dJ{v{t)) for every t G [0,T]; 

(c) the function s {q{s), f (s))jr belongs to L^{0,T), and for every 0 < t G [0,T] we have 

J{v{t)) < J{vo) + [ {q{s),f{s))jrds. 

Jo 

Remark 2.17. The explicit dependence of two constants Ci and C 2 with 

k(i)|j'' < Cl and |t'(t)|£: < C 2 for every t G [0,T], 


is given in Theorem 5.1 


Notice that in the nondissipative case, the energy inequality can not be in principle stated 
in a proper subinterval [ti,t 2 ] of [0,T] with ti > 0. This is because the measurable selection 
procedure we use to overcome lack of compactness is not in general enough to guarantee upper 
semicontinuity of the right-hand side (see Section for details). As mentioned in the Introduc¬ 
tion, a time reparametrization technique, yielding an energy equality to hold in the rescaled 
time, would allow to deal with this difficulty. 

The proofs of Theorem 2.15 and Theorem 2.16 will be given in Sectionsandrespectively. 
Before discussing them, we conclude by showing that in both cases, as a consequence of the 
stability condition, the energy equality actually holds in all the subintervals where the solution 
happens to be absolutely continuous. On the other hand, it is well-known that solutions to 
problems as those we consider here are expected to be in general discontinuous because of 
nonconvexity of the energy. 


Theorem 2.18. Let a = 0 or a = 1 and assume that the assumptions of Theorem 2.16 or of 
Theorem 2.15 are satisfied, respectively. Let v be an approximable quasistatic evolution. 

• If a = 1 and [ti,t 2 ] C [0,T], assume additionally that v is absolutely continuous in 
[ti,t 2 ]. Then 

J{v{t 2 )) +yai^iv;[ti,t 2 ]) = J{v{ti)) + [ {q{s),f{s))jrds. 


(2.14) 
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Furthermore, 

A*q{t) G dJ{v{t)) + d'ip{v{t)) for a.e. t G (ti,t 2 )- 
If a = 0, assume additionally that v is absolutely continuous in [0,i] with t>0. Then 


J{v{t)) = J{vo)+ {q{s),f{s))jrds. 
Jo 


(2.15) 


Proof. Clearly, only the ”>”-inequality in (2.14) or (2.15) has to be shown. We begin by noticing 
that, since J is locally Lipschitz by Remark 2.4, under our assumption also the map 1 1 —)• J{v{t)) 


is absolutely continuous. Let now t G [0, T] be a common differentiability point for t i—)■ f{t), 
1 1 — v{t) and 1 1 —)• J{v{t)). 

Now, for ^(t) G dJ{v{t)) we have by definition of subdifferential 

J{v{t + h)) - J{v{t)) > {^{t),v{t + h) - v{t))s . (2.16) 

If a = 0, we can take ^(t) = A*q{t) we have 

J{v{t + h)) - J{v{t)) > {q{t),f{t + h) - f{t))jr 

and differentiating 

^J(u(t)) > {q{t)J{t))jr, 

so that the conclusion follows by integration between 0 and t. 


If a = 1, we can take ^{t) = A*q{t) — C{t), with ("(t) G K* = dfj{0) in (2.16) to obtain 
J{v{t + h)) - Jiv{t)) + {C{t),v{t + h) - v{t))£ > {q{t),f{t + h) - f{t))jr. 
Differentiating we have 

d 


dt 


■Kvit)) + {C{t),v{t))£ > {q{t),f{t)) 


T ■ 


For ({t) G d'ip{0) it holds 


and therefore 


{Cit),v{t))£ < i>iv{t)) - V’(O) = V’(i'(i)) 


— J(u(f)) -h ip{v{t)) > {q{t), f{t))jr . 


(2.17) 


(2.18) 


Therefore, by integration between ti and t 2 , thanks to (2.11), we get (2.14). Now, ( 2.18| ) 
holds as an equality, and so does (2.17). Since tp{0) = 0, the inclusion C{t) G dijj{0) and 
the equality {C{t),v{t))£ = 'i/)(u(t)) together imply ({t) G dfj{v{t)). Since, by construction, 
A*q(t) = ^(t) -|- C,{t) with ^(t) G dJ{v{t)), this concludes the proof. □ 


3. Auxiliary Results 


In this section we prove some auxiliary results, that will be used to prove both Theorem 2.15 
We start by showing that, under suitable assumptions, an approximable 


and Theorem 2.16 


quasistatic evolution is automatically an evolution of critical points. 


Proposition 3.1. Suppose that ( |2.1| ), (J2), and (J3) are satisfied, and let f G iy^’^([0,T]; J^) 
and a G {0,1}. If a = 1, suppose in addition that (’^1), and hold true. Let vq ^ £ 

he a critical point of v J + aif{v — vq) in the affine space A(/(0)), and let u : [0, T] —)• <5 be 
an approximable quasistatic evolution with initial condition vo and constraint f. Then, v{t) is 
a critical point of v J + — v{t)) on the affine space A(/(t)) for every t G [0, T]. 



















LINEARLY CONSTRAINED EVOLUTIONS OF CRITICAL POINTS 


13 


Proof. Let {vsi^)k&N be as in (2.13), and let t G [0,T] be fixed. For every /c G N, let ifc G N be 
such that (to ease the notation, we do not stress the dependence of ik on t) 

^ ^ {fk P 

From the definition of approximate quasistatic evolution we have AvSf.{t) = /{ik^k)- Then, by 
continuity of / and (2.13) we obtain Av{t) = f{t). 

We thus need only to show that (5J(n(t)) + aiF*) nran(^*) / 0. By definition of constrained 
critical point, there exists qk G T', (,k £ dJ{vsi^{t)), and Ck £ d{a'ijj){vs,.{t) — V 5 ^{t — 5k)) such 
that 

A*qk =f,k + Cfc, for every fe G N. 


From (J2) it follows that J is locally bounded and therefore, by (2.13), we have 

sup J{vs^it)) < +00 . 
fceN 


On the other hand, thanks to ( |2.8[ ) and Remark 2.9, we have 

\Ck\e' < c, for every A; G N. 


Therefore, thanks to (2.1), (3.1), and (J3) 

sup \qk\T' < - sup \A*qk\£' < - ( L ( sup J(u4(t)) + 1 ) + c ) < +oo . 


fceN 


7 fcsN 


7 


fceN 


Thus, there exists q £ T such that, up to subsequences, 

lim \qk-q\T'=0- 

ft^ + OO 


From (2.13), (3.1), and (3.2) we get, by the closure property of the subdifferential, that 

A*q G dJ{v{t)) + aK* , 

as required. 


(3.1) 


(3.2) 


□ 


3.1. A constructive approach. In order to construct an approximate quasistatic evolution, 
we first introduce an auxiliary minimum problem. Let 5 G (0,1) be a fixed time step, and 
let i G N with iS < T. Set v~^ := uq, and suppose that G T is a critical point of 

u I—)• J + a'ip{v — u*“^) on the affine space A(/((i — 1)5)). If property (J2) is satisfied, we define 
the sequence (u)), j G No, by setting Vq := and 

v) := argmin{J(u) + r/|u — + aV'(u — Vq) : v £ £} for every j G N (3-3) 

Av=fiiS) 

with r] > 0 chosen such that (J2) holds. 


Remark 3.2. Note that (J2), ('I'2) and ('I'3) guarantee that mimimizers in (3.3) are unique. 
The following lemma gives some properties of the sequence (u^, j G Nq. 


Lemma 3.3. Let ( |2.1| ), (Jl), (J2), and (J3) he satisfied, let f £ 1F^’^([0,T];T"), a £ {0,1}, 
and let vq £ £ be a critical point of v J + afi{v — uq) in the affine space A(/(0)). If a = 1, 
suppose in addition that 1 (^^2) and (^'I'3 ) hold true. Let 5 £ (0,1) and let i £fi with i5 < T. 
Set v~^ := Vq, and suppose that is a critical point ofveA J{v) + afi{v — u*“^) on the affine 
spaee A(/((i — 1)5)), and let be as in (3.3). Then: 


(i) {J{v^j) + afi{v] 



is a nonincreasing converging sequence; 
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(ii) is bounded and 


lim \v) - v)_i\£ = 0] 


j—>+oo 


(3.4) 




(iii) any limit point of « critical point of the functional v i—>• J{v) + aip{v — v 

on the affine space A{f{i5)). 

Proof. For every j >2 we have Avj = Avj_i = f{i6), and therefore Vj_i is a competitor for the 
minimum problem in (3.3). Thus, 

J{vj) + aif{v) - Vq) < - Vq) - r]\vj - for every j > 2. (3.5) 

In particular, the sequence {J{vj) + a'ip{vj — Vq)) is nonincreasing. Since J + aif > 0, the 


limit 


i>2 

lim iJ{vA + aif{v) - Ug)) =: C > 0. 

j^OO ^ J / 


(3.6) 


exists and it is nonnegative, eventually showing (i). Let now M G N with M > 2. Summing up 
relation (3.5) for j = 2,..., M we obtain 
M 




Sending M 


i=2 

> oo we then have 

OO 


7 ] 


\v] - + afj{v\ - v})) - C) < OO. 


J=2 


In particular, this shows that (3.4) holds true. Note now that, by ( |3.6| ), is bounded. 

Therefore, since = \f{id)\'^ for every j > 1, the sequence {J{vj) + i® 

bounded. By (Jl), we have that is also bounded, and this concludes the proof of 


11 


Let u* G T be a limit point of Up to subsequences, we can assume that 


lim v) = u*, 

j^OO 


in 8. 


First of all, note that Tu® = f{i6). By (3.3), for every j G Nq there exists g® G if' such that 

A*q] G dJ{v]) + 2r/(u* - ?;®_i) + d{af)){v] - i;(,), 
where we used Remark 2.6 The previous relation can also be written as 

— Cj + 2??(u® — w)_i) + Cj) (3-7) 

for some ^® G dJ{vj) and Cj £ d{a'ij7){vj —Vq). Note that, since is bounded, by (J3) we 

also have that (^® is bounded. From ( |2.8[ ), is bounded. Thanks to ( |3.4[ ) and ( |2.1[ ), 

this implies that is also bounded. Thus, up to subsequences, we can assume that 


lim e® = f 

J^OO 

lim d = C* 




in 8' 


and 


lim g® = cf in T"', 


J^OO 


for some d) C* ^ and d G T"®. Passing to the limit in (3.7), thanks to (3.4) we conclude that 

A*q^ = e+c- 
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By the closure property of subdifferentials we have G dJ{v^) and C ^ d{aip)(v'' — Vq), and 
thus 

(dJ{v'‘) + d{a'ilj){v^ — Vq)) n ran(A*) 7^ 0. 

□ 

Remark 3.4. Suppose that u* and z* are two limit points of the sequence {vj)^^^^. By property 
(i) in the previous lemma and the continuity of J, even if u* 7^ z* we have 

J(u*) + — Vq) = J{z'‘) + — Vq). 

We state now a direct consequence of the previous lemma. 

Corollary 3.5. Let ( |2.1[ ), (Jl), (J2), and (J3) he satisfied, and let f G W^’^([0,T]; and 
a G {0,1}. If a = I, suppose in addition that (’^1), (^‘i) and hold true. Let S G (0,1) and 
let Vq G £ be a critical point of v J + afi{v — vq) in the affine space A(/(0)). Set v^ := vq 
and, for every i G N with i6 < T, let be defined by (3.3), and let u* be a limit point of 

{vj)j^^^. Then, the function vs '■ [—5,T] —)■ £ defined as 

vs{t) := u* for every t G [—<5, T] n [id, {i + 1)5), for every i G {—1} U No with id < T, (3.8) 

is a discrete quasistatic evolution with time step d, initial condition vq, and constraint f. 

We now prove a uniform bound for the discrete quasistatic evolution defined above. In the 
statement we write that the constant Zi depends also on ac, with c as in Remark 2.9 With 
this we mean that this additional dependence only occurs in the dissipative case a = 1. 


Proposition 3.6. Let (2.1), (Jl), (J2), and (J3) he satisfied, and let f G W^’^([0,T]; J^) 
and a G {0,1}. If a = 1, suppose in addition that (’^1), ('^‘1) and (^2>) hold true. Let 
d G (0,1) and let vq G £ be a critical point of v 1—)■ J + afi{v — vq) in the affine space 
A(/(0)). Let Vs : [—5,R] -G- £ be defined as in (3.8). Then, there exists a p ositive constant 
Zi = Zi(J,uo,7,L,ac,77, ||/||li([o,t];J-), ll/llL“([o,r];j-)), with c as in Remark 


2.9. such that 


sup \vs{t)\£ < Zi. 
<5e(o,i) 

i6[0,T] 


Proof. Let i G N with id < T he fixed, and let {v))^^^^ be the sequence defined by (3.3). By 
property (J2) and Remark 2.5 the functional „»(•) +aV'(-— Uq) is strongly convex. Therefore, 
whenever ^ G dJ^„i^{v\) + d{afi){v\ — Vq), we have 

- ^ 0 ) > + Oii’ivi - Vq) + {Lv- for every v g£. 

In particular, choosing v = Vq = v^~^ and recalling the definition of we have 

J(u*"^) > J{v\)+r]\v\ - u^ll + aTp{v\ - u(,) + - v\)£ . (3.9) 


By (3.3), v\ is the global minimizer of Jjjy^^i') + — Vq) on A{f{id)). Therefore, there 

exists r® G such that A*P G (u}) + d{afi){vl — Vq) so that, by ( |3.9[ ), 

J{v^-^) > J,y^iv{) + afiiv{ - Pq) + (AV^ Pq - vl)£, 

which gives 

Jr^,vliv\) + afiivl - v}^) < J(u®"^) + (A*r\v\-vl)}£. (3.10) 
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Step 1. We show that there exist positive constants L' = L' {L,ac,r]) and 6 = 6{L, ac,rj,'y, f) 
with the following property: for every 6 G (0, 5) and i G N with i6 <T we have 


i-G 


/ ri5 


/ riS 


|/(s)|j'ds ) J{v^) < Jiv'- ^) + — I \f{s)\jrds. 
) 7 J{i-1)5 


\ 7 

We start observing that 

AV* G + d{a'ip){v\ - vl) C + aK*. 

Therefore, thanks to condition (J3) we have 

|71V|^, <L' h.^{v\) + l), 


(3.11) 


for some positive constant L' = L'{L,ac,r]), where c is given by Remark 2.9 Thus, from (2.1) 


We also have 


\{A*r\v{-vi,)s\ = \{r\Avl-Avi,)^\ < |rV'l/(s)l^ds 


< 


< 


7/ 

7 

L' 

7 




L' 


riS 


7 / J{i-1)5 


l/(s)I.Fds 


- vh)) 


+ 


u_ 

7 


riS 




\fis)\j^ds. 


Using last inequality, (3.10) gives 


1-7 


7 Jii-1)5 


|/(s)|j-ds ) (Jrf,vl{v{) + atp{v\-vl))) <J{v^ ) + — 


L 


f piS 


7 7(i-l)<5 


|/(s)|^ds. (3.12) 


By the absolutely continuity of the integral, there exists a positive constant 6 = 6{L, ac, rj, 7 , /) 
such that 

C\G)<-5 ^ [ \f^s)\rds<\. (3.13) 

1 Jg 2 

Note now that, for every j > 2, by the minimality of u*- we have 

J{v"j) < J{v]) + v\Vj - Vj_i\£ + - u(,) < J(u*_i) + - u(,) 

< J{v\) + ailj{v\ - vl) < (vl) + ai}{v\ - vl), 

where we also took into account that {J{v'j) + a'ip{vj — vl ))is nonincreasing. Passing to the 
limit when j —)> 00 , up to subsequences, we obtain 

J{v^) < Jr,,vii'^i) + «V’(i^l - vl). 

Combining last relation with (3.12), we get (3.11). 

Step 2. We conclude. We start by proving that 

sup \J{vs{t))\ < Z 2 , 

5e(o,i) 
te[o,r] 


(3.14) 
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for some positive constant Z 2 = ^ 2 ( 7 , L', tig, ||/||ii([o,T];j'))- Note that it is not restrictive to 
assume b G (0, 5). We now set, for every i G N such that ib < T, 

Jj 

a5,i'-= and 65 ,*:=—/ |/(s)|jrds. 

7 J{i-i)S 


Notice that bs^i < - by (3.13), and that relation (3.11) gives 


^ bsi 1 

asA < -;-h 


1 — bs,i 1 — bs^i 
Iterating the previous inequality we obtain 


(^5 A — 


^ 1 ^ bgi ^ 1 


1 bs^i 1 b^^i 1 b§^{ 1 bs^i _ 1 bs^i—\ 1 bs^i—\ 


bs,i—l ^ 2 


+ 


bs,i—l 


+ 


a5A-2 


1 bs^i (1 65 ^j)(l (1 65 ^j)(l 65^4—1) 

i—l 


< ... < 


< 


a<5,o 


+ E 


bs,i—k 


(1 - bs,i){l - bs^i-i )... (1 - 65 , 1 ) ^ (1 - bs,i){l - bs^i-i )... (1 - bs,i-k) 


0-5,0 


i—l 


+ 


E'-w 


(1 - bs,i){l - bs^i-i) ... (1 - bs^i) (1 - bs,i){l - bs,i-i) ... (1 - 657 ) 


i—k 


< 


1 


<^(^ 0 )+ / |/('S)bds 

Jo 


(1 - fe<5,i)(l - ... (1 - bs^i) 

We therefore only need to find a bound for the quantity 


(1 - 65 ,i)(l - bs^i-i )... (1 - bs,i)' 


Since bs,i < -, with the elementary inequality 0 < — ln(l — x) < 2x for all 0 < x < ^ 
eventually get 

In 


(1 - 65,j)(l - bs^i-i )... (1 - 657 ) 


±Hi-bs,i)<2±bs,i<^-^ Tm 

tt ti ^ 


s)| 7 -ds 


so that (3.14) follows. We now have 


sup {J{vsit)) + \Avs{t)\jr) < sup (j{vsit)) + \\f\\looU0Tyjr))<Z3, 

<5e(o,i) <5e(0,i) ^ ^ 

te[o,r] te[o,T] 

for some positive constant Z 3 = Z^{'y,L',r], \\f\\L'^((o^T)-,j^)j II/IIl°°((o,t);J'))- Taking into account 
(Jl), last inequality implies that 

\v 5 {t)\£ < Zi, for every <5 G (0,1) and t G [0,T], 

with a constant Zi that also depends on the coercivity of the function v 1 —)• J{v) + |7lu||. □ 


4. Proof of Theorem 12.151 


We start by stating a lemma that will be used later. 
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Lemma 4.1. Let X be a Banach space, let T > 0, and let v G BV{[0,T]]X). Let C 

L°°{[0,T]] X) (1 BV{[0,T]; X) have equibounded variation. Then it exists an at most countable 
set N C [0, T] and a subsequence Skj, independent oft, with 

lim Wvs^^ft - 6kj) - vSf^Xt)\\x = 0 for every t e {0,T]\N. 

j—>-00 3 ■' 3 

Proof. We set Vfc(t) = Var(ti 5 j, ; [0,t]) and we observe that, by the assumption, is an equi¬ 
bounded sequence of monotone nondecreasing functions. By Belly’s Theorem it exists a subse¬ 
quence , independent of t, and a monotone nondecreasing function V such that 14^ (t) —)• V(t) 
for every t G [0, T], Let now N be the (at most countable) set of discontinuity points of V. For 
t G {0,T]\ N and an arbitrary e > 0, let (5o > 0 be such that 


V{t)-V{t-5o)<e. 


(4.1) 

|x), the monotonicity of 14 and 


Using the definition of pointwise variation, (2.10) (with if = 

(4.1) we get 

lim sup 11 ^ 4 , (t - 4j) - (^)IU < lim sup (^^(t) - Vk^{t - 5k^)) 

j—^oo ^ ^ 

< lim sup (14^ {t) - Vkj {t - Jo)) =V{t)-V{t-So)<e. 

j^CX) 

This proves the statement by the arbitrariness of e. 


□ 


The following a-priori estimates will be needed to prove the Theorem. 

Theorem 4.2. Let (2.1), (Jl), (J2), (J3), (’^1), (^‘2) and (^2)) he satisfied, and let f G 
1U^’^([0, T]; T"). Let 6 G (0,1) and let vq G S be a critical point of v J + ifi'a — vq) in the 
affine space A(/(0)). Let vs ■ [—J,T] -G- £ be defined as in (3.8) with a = 1. Then, there exist 
qs G L°°{[0,T];iF'), vs ■ [ 0 ,T] —>■ £, two positive constants Zi,Z^, and ri : (0,1) —)• [0,oo) and 
r 2 : ( 0 , 1 ) —>■ [ 0 ,oo) with 


lim ri(J) = lim r 2 (J) = 0 , 
( 5 —> 0 + 5 ^ 0 + 


such that 


(i) \qs{t)\x' < Z 4 and |u 5 (t )4 < Zi for every t G [0,r]; 

(ii) dist [A*qs{t) — dJ{vs{t)), K*) < 2r]ri{6), with rj > 0 as in (J2); 

(hi) sup \vs{t - 6 )-vs{t)\£ < ri{5); 
te[o,r] 


(iv) For every 0 <ti < t 2 <T 


r^2 


J{vsit 2 )) +Xar.,l,{vs-,[ti,t 2 ]) < Jivsih)) + / {qsis), f{s))jr ds + r 2 { 6 ). 


iti 


Remark 4.3. The constants Zi is given in Proposition 3.6, while Z 4 depends on the following 
quantities 

Z 4 = Z4(L,c,ry,7,'(;o, ||/||li([o,t];J-))- 

Proof. First of all, note that inequality |u 5 (t)|£- < Zi was proven in Proposition 3.6 Let now 
rj be given by (J2), let u be the constant of strong convexity of u 1 —)• J{v) + 7 |u||, and let 6 be 
given by (3.13). 


Step 1. We show that there exists a positive constant Z 4 = Z 4 (L,c, 7 , uq, ||/||l1([o,t];J')) with 
the following property. If rQ,r\ G T' are such that 

A*Pq G dJ{vi) + K* and A*r\ G dJ^y^iv{) + dX{v{ - 4), 
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then we have 


+ VWt' < ■^4- 


By (3.12) and (3.14) we have 


L' 


- Vq) < Z 2 + — ||/||l1([0,t];J-) for every 6 e (0, 5) and i E N with i6 < T, 


where Z 2 does not depend on 6 and i. Then, by (J3) and (2.1) the estimate ViIt' < Z 4 follows. 
Since by (3.14) it also holds J{vq) < Z 2 , while K* C i?(0,c) by Remark 2.9 again by (J3) and 
(2.1) the claim is proved. 

Step 2. We prove (hi). First of all, we define vs{t) = vq for every t E [0,<5), and 

vs{t) := v\ for i5 <t < {i + l)h, i E N with i 6 < T. 

Let now i E N with i5 < T. We need to show that there exists ri : [0, 6 ) —)> [0, 00 ) with ri(5) —)• 0 
as (5 —>■ 0"'' such that 

\v\-vl)\£ < ri(<5). 

By construction, we have that Vq = is a critical point of the functional v 1 —)• J(u)+i/)(u— 
on the affine space A(/((i — l)h)). Therefore, there exists Tq such that 

A*r}^ E dJ{vi^) + C + K*. 

Using the fact that cIJ(uq) = dJ^y^{vQ) we can also write 

TV* + 

Note now that, by construction, v\ minimizes the functional v 1 —>• + V’(^ “ '>^ 0 ) on the 

affine space A(/(ih)). Therefore, there exists r\ such that 

A*r\ E dJ^y^{v\) + di^ivl - u^). (4.2) 

Now, T*rQ and A*r\ are both in the subdifferential of the strongly convex functional J^^i(u) + 

ijj{v — Vq) at Vq. The constant v of strong convexity is furthermore indepedendent of Vq by 
Remark |2.5[ With this and step 1 we have 

u\v\ - u* II < (TV* - TV*, v{ - vi)e = (r* - r*, Av{ - Tu*)^ 

= (r* - r*, f{i 6 ) - fiii - l)5))r < Z^fiid) - /((i - 1)5)^ 

< Z 4 f |/(s)|j-ds. 


Then, 


\vs{t)-vs{t-6 )\‘l; = \v{-v^q\1: <— [ |/(s)|j-ds < — sup / |/(s)|j-ds. 




Setting 


i<5 \ 

ri{6) := \/— { sup [ |/(s)|j-ds) 

\i5<TJ{i-l)S J 


IZ4 


by the absolute continuity of the integral we have that ri(<5) —?• 0 as <5 —)• 0^ and we conclude. 
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Step 3. We prove (i) and (ii). We start by defining the function qs : [0, T] — )> T' . Since vq is a 
critical point of u i—>■ J + ^jJ{v — vq) in the affine space A(/(0)), there exists qq G T' such that 

A*qo£dJivo) + K*. 

We set qs{t) = go for 0 < t < <5. Moreover, we set 

qs{t) := r\ for i6 < t < (i + l)(i, i G N with iS < T. 

Then, by construction and step 1 we have \qsiiS)\j^ < Z 4 , so that (i) is satisfied. Moreover, 

A*qs{i6) G dJ^y^{v\) + 

= dJ{v\) + 2r]{v{ - + SV’K - v'o) 

caJK) + 2r?K-u^) + A*, 

so that 

dist {A*qs{t) - dJ{vs{t)),K*) < 2r]\v\ - Vq\s < 2r]ri{6). 


Step 4. We show (iv). Let 0 < ti < t 2 < T. By ( |4.2[ ) we have that for every v £ 8 

+'ipiv - vli) > Jrjy^{v\) + ipivl -vl) + {A*r\,v -v\)s, 
for every id <t < {i + l)(i, and i G N with id < T. In particular, choosing u = Uq we obtain 

J{vl) > J{v\) + i;{v\-vl) + {A*r\,vl-v\)E, 

which gives 

J{v\) + i>{v\ — Vq) < J{vq) + {A*r\,v\ — Vq)s 


nS 


= J{vs{{i - 1)5)) + 


'(*-1)5 


(®(s),/(s))j-ds. 


Thanks to (i) in Lemma 3.3 we have that for every j > 2: 


riS 


J{vj) Aipiv) - Vq) < J(vs((i - 1)5)) + 
Passing to the limit as j —>• 00 we then obtain 


'(*-1)5 


(qs(s),f(s)}jrds. 


riS 


J(vs(id)) + ip(vs(id) - vs((i - 1)5)) < J(vs((i - 1)5)) + 


'(*-1)5 


(qs(s),f(s)}jrds 


(4.3) 


Let now /, A; G N be such that Id < ti < (I -h 1)5 and A:5 < t 2 < (^ +1)5. By summing up relation 
(4.3) for i = I + 1, ■ ■ ■, k we obtain 

^ pkS 

J{vs{kd)) + '^^j;{vs{id)-vsiii-l)d)) < J{vs{ld))+ / {q5{s)Jis))jrds 


i=l 


'15 


rt2 


rti 


= J{vs{ld))+l {qs{s),fis))jrds+ I {qs{s)J{s))jrds-f {qs{s),f{s)) 
Jt\ Jl5 J kS 

<J{vsild))+ f {qsis),f{s))jrds + 2Z4Sup [ |/(s)|j'ds. 

Jti i6<T J (i—l)5 

Since vs is piecewise constant, we have 

k 

"^^Pivsiid) - vs{{i - 1)5)) = Var^(u 5 ; [ti,t 2 ]) ■ 

i=l 


T 


ds 
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With this, and since vs{k5) = vs{t 2 ) and vs{lS) = vs{ti) by construction, (iv) follows by setting 

riS 


r 2 ( 5 ) := 2^4 sup / |/(s)|jrds. 

iS<T J (i—1)5 


□ 


We can finally give the proof of Theorem |2.15 


Proof of Theorem 2.15. We divide the proof into several steps. 

Step 1. We prove that there exists a sequence and functions v € L^([0,T];T), and 

q G L“([0, r];T''), such that 

^4 —>■ n in L^([0, T]; T) and ^ Q weakly* in L°°([0, T]; T"'). 


From properties (i) and (iv) of Theorem 4.2 we have 

lk<5llLi([0T]T) + Var^(n5; [ti,t2]) < C, 

for some constant C that is independent on 6. By Belly’s theorem (and since £ has finite 
dimension), there exists a subsequence and a function v G BV{[0,T];£) such that 


and, in addition. 


n in L^([0,r];T) 


lim vs^{t) —)■ v{t) in £ for every t G [0,T] 
k^oo 


Up to extracting a further subsequence we can assume that Lemma |4.1| holds. Furthermore, 
since by (i) in Theorem 4.2 we also have 

lk<5||L°°([0,T];J'') < C, 

without any loss of generality we can assume that for the same subsequence we have 

qSk ^ Q weakly* in 
for some function q G L°°{[0,T];P''). 

Step 2. We are going to show that 

F 4 V strongly in L^((0,r);T). 

Indeed, we have 

/ \v5kis) - v{s)\£ds < / \vsi^{s) - vs^is - 6k)\£ds + \vs^^{s - 6k) - v{s)\£ ds. 

Jo Jo Jo 

By (iii) in Theorem 4.2 we have 

lim / |n 4 (s) - U 4 (s - 4)|£:ds = 0. 

k^oo In 


Concerning the second term, by Lemma 4.1, the convergence of ^4 to v in L^, and dominated 
convergence we have 

lim / |n 4 (s - 4 ) - i'(s)|£ds = 0 . 

fc->-00 In 


By (iii) in Theorem 4.2, Lemma 4.1, and the pointwise convergence of vs^. to v, we also have 

lim V 5 ^.{t) — )> v{t) in £ for every f G (0, T] \ iV , 

k—^oo 
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where N is at most countable. By (ii) in Theorem 4.2 for every t G [0,T] there exists ^s{t) G 
dJ{vs^{t)) such that 

dist {A*qs{t) < 2r}ri{5). 

In particular, since ((Z4)fcgis} is bounded in L°°([0,T];T"'), this implies that bounded 

in L°°([0,T]-,£'). Without any loss of generality, we can assume that 

^ ^ weakly* in L“([0, T]; T') 

for some ^ G L°°{[0,T]]S'). 

Step 3. We are now going to prove that 

^(t) G dJ{v{t)) for C^-a.e. t G [0,T]. 

First of all note that there exists a set A C [0, T] with >C^(A) = 0, such that for every t G [0, T] \ A 
the following properties are satisfied: 

lim - / J(v(s)) ds = J(v(t)) in M, 

h^o+ h 

lim v / (C(s),u(s))£:ds = (C(t),'(;(t))£: in M, 

h^o+ h 

1 rt-\-h 

lim - / ^(s)ds = ^(t) in T', 

h^o+ h 

lim - / |u(s)|c ds = |u(t)|^ in M, 

h^o+ h 

rt-\-h 

lim — / v(s)ds = v(t) in T. 

U^o+ h Jt 

Let now s G [0,T]. Since for every A: G N we have ^^(s) G dJ{vs^is)) = dJ^^jjg^(^s'j{vs,.{s)), 
and recalling that (s) is convex, for every w £ £ 

Jiw) > J{vSk{s)) - r]\w-vsi,is)\‘l: + (^ 4 (s),'u; - ^^(s))^. 

Let now t G [0, T] \ A. Averaging the previous inequality between t and t + h we obtain 
1 

j{w)>- [J{vs^^{s))-rjlw-vs^{s)\l + {^5 As),w-V 5^:{s))s] ds. 

Passing to the limit as fc —>■ oo we obtain (recalling that vs^. —>■ v strongly in L^) 


1 

h 

rt-^h 


rt-\-h 


- v\w - vs,,is)\1: + {^s^,is),w - V5,^{s))£] ds 


1 rt+n 1 rt+n 

^hJt + r]\w-v{s)\lds. 

Then, passing to the limit as /i —)• 0"'' we obtain 

J{w) + T]\w - u(t)|| > J{v{t)) + (|(t), w - v{t))s Mw £ £. 

The inequality above shows that ^{t) £ „( 4 )(u(t)) = dJ{v{t)). 

Step 4. We prove that 

A*q{t) £ dJ{v{t)) + K* 


for C^-a.e. t £ [0, T]. 
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Using (ii) in Theorem 4.2 we have 

pT 


dist {A*q{t) — dJ{v{t)),K*) ds 


thns giving that 


< lim inf / dist {A*(t) - dJ (t)), iiT*) ds 

Jo 

< 2r/Tliminf ri(5fc) = 0, 

fc^OO 


A*q{t) G dJ{v{t)) + K* for £^-a.e. t G [0,^]. 


Step 5. We prove the energy ineqnality. Let now 0 <ti < t 2 <T. From (iv) in Theorem |4.2| 
we have 

^(^^ 4 (^ 2 )) + Var^(r; 4 ; [ti,t 2 ]) < J{vs^{h))+ [ (qsjs), f(s)}j-ds + r 2 (Sk). 


iti 


Passing to the limit as A: —t- 00 , and nsing the fact that the total variation is lower semicontinnons 
with respect to the topology 

J{v{t 2 )) + Yai^{v, [ti,t 2 ]) < lim inf {J{v 5 ^{t 2 )) + Var^(i; 4 ; [ti,t 2 ])) 

k—^oo 


< 


liminf rj(r; 4 (ti)) + f (^ 4 ( 5 ),/(s))j-ds + r2(4) 

V Jti 


rt2 


= J{v{ti))+ {q{s),f{s))jrds. 

Jti 


□ 


5. Proof of Theorem 12.161 


We now prove Theorem 2.16 A relevant difference with the previous section is that in the 
discrete approximate energy inequality we need to choose q 5 {t) in a different way to cope later 
with the lack fo BV compactness. In particular we will need the inclusion A*qs{t) G dJ{vs{t)) 
to hold already at this level. In order to do this, we need the additional assumption (J4). 


Theorem 5.1. Let a = 0, let (2.1), (Jl), (J2), (J3), and (J4) be satisfied, and let f G 
VU^’^([0, T]; T"). Let 6 G (0,1) and let vq be a critical point of J on the affine space A(/(0)). Let 
vs : [0, T] 6 be the discrete quasistatic evolution with time step 6 , initial condition vq, and 
constraint f given by (3.8). Then, there exist a piecewise constant right-continuous function 
Q5 ■ [0, T] 


T' and positive constants Ci, C 2 and C 3 , independent of 6 , such that 

(i) A*qs{t) G dJ{vs{t)) for every t G [0,r]; 

(ii) \qs{t)\T' < Cl and \vs{t)\£ < C 2 for every t G [0,r]; 

(Hi) for every ti < t 2 G [0,T] 

J{vs{t 2 )) < J{vs{ti)) + f 

Jtx 


{qs{s),f{s))T<^s + 


Remark 5.2. More precisely, as it appears by a careful reading of the proof of Theorem 5.1 
we have 


Cl — C'i( J, 7, L, r?, ||/||l 2 ((o,T);J')) a-Ild Cs — C3{J,'J,L,r],Cj^ri, II/IIi,2((0,T);J'))- 
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The boundedness of v{t) has been already proved in Proposition 3.6 and one can take for C 2 
the constant Zi provided there. The dependence on J has again to be unde rstoo d in the sense 
of the coercivity of the function u 1 —)• J{v) + see the proof of Theorem 


4.2 


Proof. Since vs is a discrete qnasistatic evolution, for every i G N with i5 < T there exists 
g* G F' such that A*q'‘ G dJ{vs{i5)). Then, if we define qs : [0,T] —)■ F' as 

qs{t) := q'‘ for all t G [0, T] H [i5, (i + 1)(5), for all f G Nq with i 6 < T, 

property (i) is satisfied. Since |u 5 (f)|j- is bonnded by Proposition 3.6 with (J3) and (2.1) we 
obtain (ii) We now divide the remaining part of the proof into two steps. 

Step 1. We show that there exists a constant M, depending only on rj and such that 

JivsiiS)) < J{v 5 {ii-l) 6 ))+ [ {q\s)J{s))jrds + M 6 [ \f{s)\‘jrds, (5.1) 
for every i G N with i5 < T. 

To this aim, let i G N with i 6 <T he fixed, and let be the seqnence defined by (3.3). 

By property (J2) and Remark 2.5 the functional strongly convex. Therefore, whenever 

i G we have 

Jr^,vi{v) > Jr,,vi{v\) + {i,v - v\)s for every u G T. 

In particular, choosing v = Vq and recalling the definition of we have 

|2 


> J{v\) + r]\vl - 


for every ^ G (5.2) 

By (3.3), v\ is the global minimizer of on A{f{i5)). Therefore, there exists r* G F' such 
that A*P G so that, by ( |5.2[ ), 

> J (uj^) + r]\v\ — UqII + {A*r\vQ — v\)£. 

Therefore, by the absolute continuity of /, and recalling that qs is constant in the interval 
[(i — 1)(5, i(5), we have 

rj\v\ -vi\l + J{v{) - J(u*) < {A*r\ vl - vi)s 
= {A*P - A*qs{{i - 1)<5), - Vq)£ + {A*qs{{i - 1)(5), - vlf)£ 

= y - qsiii - 1)6), Avl - Avl)jr + {q\{i - 1)5), Av\ - Avl)jr (5.3) 

= 'f'' - qsiii - 'i-)6),Avl - Avl)jr + f {q\ii - 1)(5), /(s)) ds, 

where we used the fact that Av\ = fii6) and Avq = /((z — 1)5). Observe now that, by definition 

of qs, we have A*qsiii — 1)5) G dJivsHi — 1)5)). Thus, recalling that vsHi — 1)5) = = Vq, 

we obtain 

A*qsiii - 1)5) G aj(u^) = dJ^„^ivl). 

Thus, recalling that A*P G dJ.^y^iv{), by property (J4) we obtain 

(P - qsiii - 1)5), Av\ - AvDjt < Cj^nWi - vl\£\Av\ - Av^fjr 
which, together with ( |5.3[ ), gives 

V\v\ - vl\l + Jiv\) - Jivl) < Cj^nWi - u(||£-|/(z 5) - /((z - 1)5)17- + f {q\s), f is))jr ds. 

J (i—1)5 
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Using Young’s and Holder’s inequality, we get 




JK) < J{vl) + M\f{i5) - fiii - l)5)|j. + / {q\s), fis))^ds 

2 


= + M 


piS 

J fz—] 


< J{vi)) + M 6 


fis)ds 

iS piS 

2 


+ [ {q^{s), f{s))jrds 

J- 

riS 

S( 


(5.4) 


' |/(s)|_^ds+ / {q {s)J{s))jrds, 


for a suitable constant M > 0 (depending only on r/ and where we also used the fact that 

/ G VU^’^([0,T]; J"). Recalling that J{vj) < J{v\) for all j > 2, by (3.8) and property (ii) in 
Lemma [3.3t we have 

J{vs{iS)) = J(f*) = lim J{vj) < J{v\). 


J -!-00 


Taking into account last inequality, and recalling that Uq = n* ^ = vs{{i — 1)5), relation (5.4) 
gives (5.1). 

Step 2. We show (hi). Let G N be such that 16 < ti < {I + 1)5 and k6 < t 2 < {k + 1)5. 
Recall that vs{kS) = vs{t 2 ) and vs{l6) = vs{ti) by construction. By summing up relation (113 
for z = ^ + 1,..., A:, using Holder’s inequality and taking into account that 5 G (0,1) we have 

pkS 

Jiv5it2)) < J{vsih)) + I (^5(s),/(s))j-ds + M5||/|||2((0,T);.F) 

r^kS 


ns 


pt2 _ pti _ pKO 

J{v5iti))+ {q5{s),f{s))jrds+ {qs{s)J{s))jrds- (/(s),/(s)) ds + M5||/||^2 ((o,t);J-) 
J t\ J IS J t2 

< J{v5{ti)) + f {q^(s), f (s))jr ds + C'i'v/5||/||L2((o^'r).j-) + lH\/5||/||^2((o^r);j') 

J tl 

= J{v5{ti)) + [ {q^{s), f{s))Tds + Vs (^C'i||/||L2((o,r);j') + -^ll/lli2((o,r);j^)) > 


which gives (iii). 


□ 


Before giving the proof of Theorem 2.16, we need the following result (see |12l Lemma 3.6]). 

Lemma 5.3. Let X he a compact metric space. Let p : [0,r] —>■ M, pfc : [0, T] —>■ M and 
/fc : [0, r] —)• Y be measurable functions, for every fc G N. For every t G [0, T] let us set 


X{t) := {x G Y : d/cj —5- +oo such that x = lim fkjV) and p{f) = lim pfc^(t)}. 


J—>+CXD 




Then 


• X{t) is closed for all t G [0,T]; 

• for every open set U C X the set {t G [0, T] : X{t) nU 0} is measurable. 

We conclude this section with the proof of Theorem |2.16[ 

Proof of Theorem 2.1()[ We divide the proof into several steps. 

Step 1: proof of (a) and (b). 

For every 5 G (0,1), let qs : [0, T] —> X' and vs ■ [0, T] —^ T be given by Theorem 5.1 Let 
A C [0,r] be such that T^(A) = 0 and f{t) is well defined for every t G [0, T] \ A. We fix a 
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sequence such that Sk —t 0"'' and define 


9k{t) := 


for every t G [0, T] \ A, 
for every t G A, 


and 


6{t) := limsup0fc(f) for every t G [0,r]. 
k^oo 


By definition of 6, for every t G [0, T] we can extract a subsequence (possibly depending 

on t) such that 

9{t) = lim 6k.{t) for every t G [0,r]. 

j^oo ^ 


By (ii) of Theorem 5.1 we have 

\qs{t)\T' < Cl, \vs{t)\£ < C 2 , for every (5 e (0,1) and t G [0,r]. 

Thus, for every t G [0, T] we can extract a further subsequence (not relabelled) such that 
lim (t) = v{t) in S and lim (t) = q{t) in T"', 

j—^oo j j—^oo J 

for some v{t) G £ and q{t) G T' with \q{t)\jri < Ci and |u(t)|£: < C 2 - Let us now show 
that, for every t G [0, T], we can choose the subsequence in such a way that the maps 

q : [0, T] —)• T' and u : [0, T] —)• T are measurable. 

Let us denote by (^ 02 ) closed ball of T" {£) with center at the origin and radius Ci 
{ 02 )- Applying Lemma 5.3 with X = x B^^, fk = (qs^j'^s^) pk = Ok, we have that 

• X{t) is closed for all t G [0,T], 

• for every open set U X X the set {t G [0, T] : X{t) n C/ 7^ 0} is measurable, 
where the set X{t) is given by 

X{t) := {{q{t),v{t)) G Bq^ x Bq^ : 3 kj ^ +00 such that 

{q{t),v{t)) = .lim (g4,(t),U4,(t)) and e{t) = lim OkAt)}. 


j—>-+oo 


J^OO 


Thanks to [8l Theorem III.6], for every t G [0, T] we can select {q{t),v{t)) G B^'^ x B^^ such that 
t —)■ {q{t),v{t)) is measurable. Thus, (a) is proven. Finally, by repeating the arguments used in 
the proof of Proposition 3.1 we obtain (b). 

Step 2: Proof of (c). Observe that, for every t G [0, T] \ A, 

9{t) = lim sup ^^(t) = lim 9k.{t) = lim {qs^^Xt), f{t))jr = {q{t), f (t)) jr. 
k^oo ^ 

Let us now show that 9 G L^{0,T). Since 9 is the lim sup of measurable functions, we deduce 
that it is measurable. Moreover, we have 

r m\dt = r \{q{t), f{t))r\dt < Cl r dt < ci%/r ii/iii2((o,T);.F)- 

Jo Jo Jo 


In order to get the energy inequality, recall that by (hi) of Theorem |5.1| (for t 2 = t and ti = 0) 
we have, for every j G N, 


^(^'4,(^)) < ^(^^o) + / (94.('®)j/(s)).^ds + C3 5 


1/2 






LINEARLY CONSTRAINED EVOLUTIONS OF CRITICAL POINTS 


27 


Taking the limsup in j of the previous expression, using Fatou’s Lemma 


= lim J{vsi,.{t)) < J(no) +limsup / (qs^Xs), f{s)) ds 


J^OO 


j^oc Jo 


rt 


< J(no) + limsup / {qs^{s), f{s))jr ds < J{vo) + limsup /(s))j-ds 


k^oo Jo 

= J{vo)+ [ {q{s),f{s))jrds, 
Jo 


/ 0 fc^oo 


so that (c) follows. 


□ 


6 . A DISCRETE VERSION OF FRACTURE EVOLUTION FOR COHESIVE ZONE MODELS 

In the remaining part of the paper, we show how to apply Theorem |2.16 to cohesive fracture 
evolution. We start this section by recalling the model introdnced in [ 6 ], where a critical points 
evolution is obtained by following the scheme described in the first part of the Introduction. We 
conclude the section performing a finite dimensional discretization. In the next section we will 
then show how it is possible to pass to the limit, thus obtaining a different proof of the existence 
result in [ 6 ]. 


6.1. A previous model for the time evolution of cohesive fractures. Let C be 

a bounded open set with Lipschitz boundary, with d G {1,2}. We assume that the reference 
configuration is the infinite cylinder x M C and that the displacement U : 11 x M —)> 

has the special form 

U{xi,... ,Xd,Xd+i) = (0 ,...,0,u(xi,...,xd)), 

where u : H —?• M. This situation is referred to in the literature as generalized antiplanar shear. 

We assume that the crack path in the reference configuration is contained in T n H, where 
r C is a Lipschitz closed set such that 0 < n H) < oo and H \ T = H"*" U 11~, where 

H"'' and 11“ are disjoint open connected sets with Lipschitz boundary. 

We will study the energy of a finite portion of the cylinder, obtained by intersection with two 
horizontal hyperplanes separated by a unit distance. Given a time interval [0, T], with T > 0, we 
assume that the evolution is driven by a time-dependent displacement oj G iL^([ 0 , T]; Ll^(ll)), 
imposed on a fixed portion of (9H. We make the assumption that is well separated 
from T, and that n 911^) > 0 . 

In the framework of linearized elasticity, the stored elastic energy associated with a displace¬ 
ment u G \ r) is given by 

W{u)-.= \l \VuXdx. (6.1) 

2 Jq\t 

The crack in the reference configuration can be identified with the set 

:= {x G T : u'^{x) / u“(x)}, 

where denotes the trace on T of the restriction of u to H^. In order to take into account 
the cohesive forces acting between the lips of the crack, according to Barenblatt’s model [5] we 
consider a fracture energy of the following form: 
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where [rt] := — u~ is the jump of u across F, and the energy density g : [ 0 , oo) —)• [ 0 , oo) is 

a , nondecreasing, bounded, concave function with < 7 ( 0 ) = 0 and a := g{0'^) G (0,+ 00 ). We 
will consider here the special case (see also [HI Section 9]) 


9{s) = < 


R 

2 


if 0 < s < R, 


if s > R, 


( 6 . 2 ) 


where the parameter R > 0 represents the range of the cohesive interactions between the lips of 
the crack. The parameter k is the stiffness material constant and, for sake of simplicity, it will 
be taken equal to 1 through all the theoretical analysis. Summarizing, the energy functional 
E : ^ [0, 00 ) is given by 


E{u) = l [ |Vu|^da:+ [ g{\[u]\)d'H'^ ^ 

^ Jn\r Jr 

Let t G [0,T], and let u{t) be a minimizer for the problem 

min{iii(u) : v G H^{Q \ F), u = Lj{t) on Sun}. 
Then [HI Proposition 3.1], 

\7u{t)-V'tpdx+ [ ([V^] 5 '(|[M(t)]|)sign([u(t)])lj + |[V^]|lj. 

' Jr ^ 


(6.3) 


/0\F 




(6.4) 


for all '0 G \ F) with -0 = 0 on 5^)0. One can see | 6 [ Proposition 3.2] that u{t) satisfies 

(6.4) if and only if it is a weak solution of 


Au{t) = 0 
u{t) = uj{t) 
dyu{t) = 0 
duU^{t) = dyU~{t) 
\dyu{t)\ < 1 


in Q \ F, 
on 

on \ SdO, 
on F, 

r \ Ju(t) 1 


(6.5) 


dyu{t) = 5 '(l[u(t)]l)sign([u(t)]) on 


F{u,t) = 


where sign(-) denotes the sign function. Any function u{t) satisfying (6.4) or (6.5) will be referred 
to as critical point of E at time t. 

In | 6 ], a critical points evolution is obtained by following the general ideas given in the 
Introduction, by setting Y = and 

E{u) for u ^ H^{U,\V),u = w{t) on 
+00 otherwise in L^(fl). 

Given a critical point uq of E at time 0 and £ G (0,1), the author shows [HI Theorem 4.1] the 
existence of a function : [ 0 ,r] —)■ satisfying m^( 0 ) = uq and ( 1 . 2 ), which in this case 

reads as 

Au^{t) = eif {t) in n \ F, 

u^{t) = ijj{t) on 

dtjU^{t) = 0 on cin \ 

duU%t)\n- = duU\^+{t) on F, 

\d^u%t)\<l onF\J„^(j), 

= ff'(l[u^(t)]l)sign([u^(t)]) on 


( 6 . 6 ) 
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Under the additional assumption g G (which is fulfilled by the function in (6.2)), uniqueness 


for the above problem also holds true [6l Theorem 4.2], Finally [6l Theorem 4.4], there exists 
a bounded measurable function u : [0,r] —>■ H^{^} \ T) with u(0) = uq such that the following 
properties are satisfied: 


approximability: for every t G [0, T] there exists a sequence 

u{t) weakly in H^{Q \ T); 


0"'' such that 


stationarity: for a.e. t G [0, T] the function u{t) is a critical point of E at time t; 
energy inequality, for every f G [0,T] 

E{u{t)) < E{u{Q)) + f [ Vm(s) • Va;(s) dxds. 

Jo Jq\v 


(6.7) 


6.2. Discrete Setting. In view of the applications of the results of this paper to the model 


just introduced, we need a finite dimensional version of the energy functional E in (6.3). In 


order to focus on the main ideas of our approach, we keep the formulation as clear as possible, 
considering a very simple geometry. 

Let d = 2 and £ > 0 be fixed, and define 

II:=(0,2t)^ Q- := (0,7) X (0,2^), Q+:= (7, 27) x (0, 27), T := {7} x [0, 27]. 

We will study a fracture evolution where the deformation is imposed on the set 

:= ({0} X [0, 27]) U ({27} x [0, 27]), 

see Figure 


(0,27) r (27,27) 



Figure 3. Geometry of the problem. 

For every h > 0, we assume we are given a triangulation Th of the set 12 \ F, and we dehne £h 
as the finite dimensional space of continuous functions that are affine on each triangle belonging 
to Th- More precisely, we set 

£h '■= {u G C{Q \ F) n \ F) : Vu = const, a.e. on T, for every T G Th} , 
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and we define as the set of functions of £h that do not jump across F: 

■.= £hr\H\^l). 

We endow £h with the induced norm of H^{£l \ F) 

:= / tt^dx+ / |Vrt|^dx u & £h, 

Jn Jn\r 

With a slight abuse of notation, we will denote with (•, ■)£^ both the duality pairing between £h 
and £h (whenever the two are not identified), as well as the scalar product in £h. In particular, 
throughout this section, we convene that the equality 


C = v 

where ^ G £h and v G £h, is meant in sense of the Riesz isometry. We will denote the restriction 
of the energy functionals E and W to the space £h by 


Eh ■■=E 


Wh := W le, . 


We denote by Ah the operator which associates to every function of £h its trace on doEf-, and 
we set Fh '■= Ah{£h)- Note that Fh is closed, since £h is finite dimensional. Therefore, Fh 
endowed with the induced scalar product is a Hilbert subspace of We also 

notice that, since the operators Ah are the restrictions to £h of the surjective trace operator 
A : —)• by the open mapping theorem we can assume that ( 2 . 1 ) is satisfied 

for all h with a constant 7 independent of h. In the sequel this will be used tacitely. 

We finally observe that, applying m Lemma 4.1.3] to our setting, we obtain the following 
version of Poincare inequality: 


Ik - «d||l 2 (q) < C'||Vu|k 2 (Q\r) for every u G £h, 
where the constant C depends on H and and 



The previous inequality in particular implies that 

IkllLRo) <UD + C'||Vu|k2(n\r) < C {\Ahu\jr^ + ||Vu||i 2 (f^\r)^ for every u G £h, (6.8) 

where with C we denote different constants, all depending on £l and We conclude this 

subsection with an important remark that will be used later. 


Remark 6.1. Let v G £h, w G £h^^, and let ^ G dEh{v). Then, from the definition of 
subdifferential and direct computation, one can check that the action of ^ on tc coincides with 
the action of the Frechet differential dWh{v) on w. In formulas: 

= {dWh{v),w)E^. 


We show now that the functional Eh satisfies the assumptions of Theorem 2.16 


6.3. Assumptions (Jl)—(J4) are satisfied by Eh and Ah. First of all, we start by observing 
that condition (Jl) is satisfied, by using standard arguments of calculus of variations. 


Proposition 6.2. The functional 

£h 3 V I —> Eh{v) + 

is coercive. 
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Proof. Let C > 0 be fixed, and let C Sh be a sequence such that 

Eh{vk) + < C- 


Then, recalling the expression of Eh and thanks to Poincare inequality (6.8), we have 

\vk\sf, < C, 

for some new constant, still denoted by C, depending on Then, there exists a subse¬ 

quence {vkj)j^^ and a function v ^ Eh such that 

Vkj V weakly in Eh- 
Since Eh is finite dimensional, this implies that 


Vkj V 


in Eh 


and this concludes the proof. 

We now show that condition (J2) is satished. 

Proposition 6.3. There exists rj > 0 such that the function 

— > Ehiu) + ri\u\l^ 

is strictly convex. 

Proof. We divide the proof into several steps. 

Step 1. We show that there exists // > 0 such that the function 


□ 


given by 


:=5(|s|) +h■s^ 


s G 


(6.9) 


is strictly convex. To this aim, we need to hnd p such that the second distributional derivative 
p'h of Pfj. is a positive Radon measure. Recalling the definition of g, we have 


»i + u-57; i»‘ 
p4») = {„ ^ ' 

Tl 9 

The distributional derivative of is given by 

Pu(s) = < 


if 0 < |s| < R, 
if |s| > R, 


—1 -|- ( 2p — — ) s if — R < s < 0, 


2fis 


if 0 < s < R, 
if |s| > R. 


Note that p'^ G We can then calculate the second distributional derivative p” of p^, 

which is the Radon measure in M given by 

P'h = ^ L(-R,R)+2h-^^ L(-oo,-R)U(R,oo)+2(5o, 

where <5o represents the Dirac measure concentrated at the origin. Note that 

p'li{B) > ^2p — C^{B), for every Borel set R C M. 
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Thus, if we choose /i such that 

/^ > ^, (6.10) 

p"^ is a positive Radon measure on M, and is convex. 

Step 2. We show that the functional ■ £h [0, oo) given by 

Eh{u) := Eh{u) + p j^\[u\\^ , 

is convex. By the previous step, the function ^ [0, oo) defined as 

r^li[u]) := J^g{\[u]\) dn^ + p dn^ 

is convex, where Qh is the subset of T^(r) given by 

Gh := {[u] : u G £h}. 


Note now that 


Eh{u) = Wh{u) + rf,{\[u]\), 


where W was defined in (6.1). From the fact that Wh ■ £h [0, oo) is convex, we then obtain 
that also E^ : £h —)• [0, oo) is convex. 

Step 3: conclusion. By [6l Lemma 5.3], there exists a constant C > 0 such that 

j^\[ufdn^<C\uW. ( 6 . 11 ) 

Taking p > pC we have 

Eh{u) + T/|u||^ = Eh{u) + p J^\[u]f dx 2 +p (^C |m||^ - ^ IMP dx 2 ^ +(r? - pC)\u\l^. 


Eh{u) 


Eh{u) 


We have already proven that E^ is convex. On the other hand. Eh is a quadratic form which is 
positive semidefinite by (6.11), and thus is convex. Since the remaining term (p — pC)\u\£^ is 
strictly convex, this concludes the proof of (J2). □ 


Before passing to the proof of (J3) we need some preliminary results. First, we make a few 
remarks on the regularity of the elastic part and on the crack part of the energy functional. 

Remark 6.4. Note that Wh G (^^^’^(T/i). In particular, dWh{-) : £h ^ £h is a single-valued 
Lipschitz function with Lipschitz constant 1. Indeed, we have 


{dWh{w),v)£^ = / Vw-Vvdx for every tt), u G T/i. (6-12) 

Jn\r 

Then, for every wi,W 2 ,v G £h 

\dWh{wi) — dWh{w 2 )\£/ = sup < / (Vwi — Vw 2 ) ■ Vvdx, v £ £h with = 1 r 

[Jn\r J 

< sup{||Vu;i - Vrn211^2IIVu 11^2 v £ £h with = l} 

< sup{|u;i - W 2 \£jv\£^, V £ £h with = 1} 







LINEARLY CONSTRAINED EVOLUTIONS OF CRITICAL POINTS 


33 


Since dWh{^) = 0, this implies 

\dW{w)\^i < \w\£^ 1 for every w G 8h- 
Remark 6.5. From the previous remark, it also follows that 

\dWh{wi) - dWh{w 2 )\£i^ < \\Vwi - Vir2||L2(D\r) for every wi,W 2 G £h- 
Remark 6.6. The functional Gh ■ £h ^ [0; oo) defined as 

Ghiv) := J^g{\[v]\)dn^. 

is globally Lipschitz continuous. Indeed, for every vi,V 2 G Sh we have 

|G',(ui) - G,{v 2 )\ < ^ |5(INI) - 5(INI)| ^ ||[i;i]| - \[v 2 ]\\dn^ 

< j \[vi] - [V 2 ]\ dn^ = j \[vi-V 2 \\d'H^ 

< IKui - V2]\\LHr) < C |i;i - V 2 \£„ 


(6.13) 


where G is given by (6.11), and we used the fact that ||5 ^||l°o([o,oo)) = 1- 
Next proposition shows condition (J3). 

Proposition 6.7. Efi satisfies condition (J3). 

Proof. Note that 

Eh{u) = Whiu) + Gh{u), 


where Gh ■ £h —^ [O^oo) is defined in Remark |6.6[ Let now u G T. By (2.6), every ^ G dJ{v) can 
be written as 

f, = f,i+ f,2, 

where G dWh{v), and .^2 £ dGh{v). By Remark |6.5| we have 

< ||Vu||2,2(Q\r) < 1 + ll^r;||^2(Q\r) — + ^h{v)). 

On the other hand, thanks to Remark |6.6| Gh is globally Lipschitz continuous with Lipschitz 
constant G Therefore, by (2.5) 

\f2y^<c{H\r)f\ 

lei^' < lak;. + 161^' < 2Ehiv) + 2 + O (Ri(r))'/' . 

□ 


Thus, 


Next lemma will be used to prove (J4), and gives a bound on the norm of a regular critical 
point, in terms of its trace on do^i- 

Lemma 6.8. Let w G £h^^, let f G Eh be such that AhW = f, and suppose dWh{w) G ran(^*). 
Then, there exists a positive constant G = C{Ll,^) such that 

\w\sh < C'l/b^- 


(6.14) 
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Proof. By ( |6.8[ ) and ( |6.12[ ) 

IMsh ^ < C [\Ahw\^j:^ + ||VR;||^2(Q\r)) (6.15) 

= C[\f\'jFt, +{dWh{w),w)E^], (6.16) 

where C denotes different constants, depending only on fl and Sufi. Let now g G be such 
that A\q = dWh{w). Then, 

\dWh{w)\^, = \Alq\^^, > 7|g|j-^. 


Thus, taking into account (6.13) we have 


{dWh{w),w)s^ = {Alq,w)s^ = {q,Ahw)jr^ = {qj)^^ < |g|j-J/| 


J^h. 


< i \dWhiw)\£,^ i \w\g^ 


Using (6.16), we obtain 


^ 


Ik + “ 


From the previous relation, the thesis follows using Young inequality. 
We can finally prove (J4). 

Proposition 6.9. Let ui,U 2 ,u G £h o.nd qi,q 2 , fi, f 2 £ J^h be such that 

AjiVi — fi, ^ d{,Phfr),v{pi)i i — 1 ) 2 , 


□ 


where rj is given by Proposition 6.3. Then, there exists C > 0, depending only on 11 and 
such that 

{Q1 — Q2,AhVi — AhV2)Th ^ C!\vi — V2\£f}AhVi — AhV2\Th- 
Proof. Let w be the unique solution of the following minimization problem: 


w = argmin{lU/i(u) : Ahiv) = fi - f 2 } ■ 


vGSr 


(6.17) 


By Remark 6.1 we have 


{f.i,'w)£^ = {dWh{vi), w)£f^ i = l, 2 . 

Now, by definition of and Remark |2.6| there exist G dEh{vi), with 1 = 1, 2 such that 

- 2T/(uj - u) = z = l,2. 

Subtracting term by term we obtain 

^h{<P - 92) - 2 ry(ui - V2) =ii- f 2 - 

Thus, thanks to Remark |6.4| 

(gi - q2,AhVi - AhV2)rh = {<li “ 92 , h - f2)Ty, = (91 - q2,Ah{w))jr,^ 

= (^h( 9 i - 92), = (6 - i2,w)£^ + 2 r]{vi - V2, w)£^ 

= {dWh{vi) - dWh{v2),w)£f^ + 277(ui - V2,w)£,^ 

< (l + 2r/)|u;|£-Jui - U 2 |£-^ < 0(1 + 2r/)|/i -/ 2 |j-Jui - U 2 |£-;,, 


where we also used the fact that w satisfies the assumptions of Lemma 6.8 with f = fi — f 2 - D 

We conclude this section with a remark, which clarifies why we didn’t directly prove the main 
theorem of the paper in the infinite dimensional setting. 
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Remark 6.10. As already explained in the Introduction, the proof of the main result of the 
paper (Theorem 2.16) is given only in a finite dimensional setting. This is due to the fact that, 
in the infinite dimensional case, the subdifferential is in general not closed with respect to the 
weak convergence in the domain of the energy. Such a difficulty could be overcome by requiring 
that the energy functional has compact sublevels, an assumption which is quite common in 
literature. In the cohesive fracture model case, this would amount to choosing as domain 

of the energy , and considering the Dirichlet linear constraint as being encoded by an unbounded 
densely defined surjective linear operator A : L^{VL) —)• This does not affect neither 

(2.1) nor hypothesis (Jl), which are still satisfied with minor modifications. Even condition (J2) 
can be proved with a little bit more of effort. On the other hand, the key conditions (J3) and 
(J4), which we need in order to cope with the cohesive fracture energy in our model, would no 
longer hold true. This motivates our choice of first dealing with a finite dimensional setting, and 
then extend the results with a problem-specific technique. 


7. Recovering an approximable quasistatic evolution 


We show now that the existence of a quasistatic evolution for the functional E, in the sense 


of [6| and Section 6.1 can be recovered from a discrete quasistatic evolution for when the 


parameter h controlling the mesh size tends to 0. Before stating the main theorem of the section 
we need some notation. We set 


D :=■{ h > 0 : h = — for some iV G N 


with £ as in 6.2 Let oj G W^’^([0, T], R^(R)), and let uq be a contrained critical point of E at 
time 0, under the constraint uq = w(0) on By [M], there exists a sequence {u:h)h^£) C 

such that Uh G Sh) for every h £ D and 

CO (7.1) 

For every t G [0, T], we define fhit) ■= AhUhit), we will denote by f{t) the trace of Lo{t) on 
Again by [3l|, there exists a sequence {uQ^h)f^^£, C \ T) such that 

uo^h G £h with AhUQ^h = /h(0) for every h G T>, 

and 

iniL^(R\r). (7.2) 


/i-s>0+ 

^ Uq 


Remark 7.1. Let h £ D he fixed. Since in general no,/i is not a critical point of E^ at time 
0 , it is not possible to consider an approximable quasistatic evolution with initial condition 


uo,/i and constraint fh- We can, however, modify Dehnition 2.14 in such a way that no critical 


point condition is required at the initial time 0, as clarified below. For our purposes, this is 


still sufficient. Indeed, as stated in Theorem 7.5 below, the critical point condition at t = 0 is 
recovered when passing to the limit /i —)• 0. 


Definition 7.2. Let h £ D, 5 £ (0,1), and let uo,/i fh be given above. A discrete quasistatie 
evolution with time step 5, initial condition uo,/ij constraint fh is a right-continuous function 
us^h ■ [0, T] —)• Sh such that 

• 1^5,/i(0) = Uo,h] 

• us is constant in [0, T] n [i<5, {i + 1)(5) for all i G No with id < T; 

• us{i5) is a critical point of Eh on the affine space A{fh{iS)) for every i G N with id < T. 
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Remark 7.3. Note that, in the definitions above, we do not require us^h{0) to be a critical point 
of Eh at t = 0 . 


Definition 7.4. Let h € D, and let uo,/i and fh be given above. We say that a measurable 
function Uh '■ [0, T] —)• £h is an approximable quasistatic evolution with initial condition uo,h and 
constraint fh, if for every t G [0, T] there exists a sequence 6k —>• 0^ (possibly depending on t) 
and a sequence {us,.^h)k£^ of discrete quasistatic evolutions with time step 6k, initial condition 
uo^h, and constraint fh, such that 

lim \us^^h{t) - Uhit)\£h = 

/c^+oo 


With this choice, an approximable quasistatic evolution still satisfies Uh{0) = uo^h- On the 
other hand (see Theorem 7.6 below) in this case we will simply require that the stationarity 
condition for Uh holds for every t G (0,T], while stationarity at 0 will be recovered in the limit 
passage /i —)■ 0. 

Our goal is now proving the following version of [6l Theorem 4.4], which is the main result of 
this section. 


Theorem 7.5. Let oo G and let uq be a eritical point of E at time 0 with 

uq = a;(0) on OdLI- For every h £ D, let ojh and uo^h be defined as above. Then, there exists 
a bounded measurable function u : [0,T] —)• \ T) with u(0) = uq such that the following 

properties are satisfied: 


(a) approximability: for every t G [0,r] there exists a sequence hj — )• O'*" such that 

Uhjit) u{t) weakly in H^{Ll \ T) 

where, for every j G N, Uh^ is an approximable quasistatic evolution of Eh^ with initial 
condition o,nd constraint fh,,, see Definition 1.2: 

(b) stationarity: for every t G [0, T] the function u{t) is a critical point of E at time t under 
the constraint u{t) = Lo{t) on 

(c) energy inequality: for every t G [0,T] 

E{u{t)) < E{u{0)) + f f Vtt(s) • Va;(s) dxds. 

Jo Jn\r 


Before proving Theorem 7.5 we need the following result, which is obtained by applying 
Theorem |2.16| to Eh. 


Theorem 7.6. Let h £ D be fixed, and let and fh given above. Then, there exists a 
measurable bounded mapping Uh'. [0,T] — )■ £h sueh that 

(a’) Uh{-) is an approximable quasistatic evolution with initial condition uo^h o.nd constraint fh; 

(b’) stationarity: for every t G (0,T] we have 



Vuh{t) ■ Vip dx + 


[i’]9'{\[uhit)]\)sign{[uh{t)])lj^^^^^ + dn'^ ^ > 0, (7.3) 


for all € £h with = 0 on doD. 
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(c’) energy inequality: The function s 


Vuh{s) ■ Vujh{s) dx belongs to L^{0,T) and 


Eh{uh{f)) < E{uo^h) + 


Jn\T 

Vuh{s) ■ VCjh{s) dxds 


for every t G [0, T], 


/o Jn\T 

(d’) Uniform bound: There exists a constant C 2 , independent of h, such that 

\\uh{t)\\£^ < C 2 for every t G [0,r]. (7.4) 

Remark 7.7. As already observed, property (a’) has to be intended in the sense of Definition|7.4[ 


Proof. As proven in the previous subsection, jg and assumptions (J1)-(J4) are satisfied. We 
now need to check that the proof of Theorem 2.16 can be adapted to the present situation, where 
Definition 7.2 and Dehnition 7.4 substitute Definition 2.13 and Dehnition 2.14, respectively. 

Step 1 : We construct a discrete quasistatic evolution. Let 6 G (0,1) be a hxed time step, and 
let i G N with id < T. We set := while, for z > 2, we suppose that G is a 
critical point of E^ on the affine space A{fh{{i — 1)5)). Then, analogously to (3.3), we dehne 
the sequence (u)jj) by setting := and 


D’GNo 

:= argmin {Eh{v) + r]\v - < .iH : v £ Sh} 
AhV=fh{iS) 


for every j G N. 


Let now z G N with id < T. One can check that Lemma 3.3 and Corollary |3 . 5| still hold true. In 
particular, for every z G N with id <T we can find a function u\ G Sh with Ahu\ = fh{id) such 
that, up to subsequences, 

lim u\ = u\ in Eh- 

j^OO 

Moreover, the function Uh ,5 '■ [0, T] —)• Eh defined as 

Uh,5{t) ■= for every t G [0, T] 0 [z5, (z + 1)5), for every z G No with id < T, 

is a discrete quasistatic evolution with time step 5, initial condition {uo)h, and constraint fh, in 
the sense of Definition 17.21 


At this point, we need to check that Theorem 5.1 can still be proven. In particular, we want 
to define a function qh^s '■ [0, T] —>• Fh such that an approximate energy inequality (as (iii) 
of Theorem 5.1) holds true. The main problem consists in defining qh^s in the interval [0,5). 
Indeed, since zi/i,5(0) is not a critical point, we do not have a natural choice available. We then 
modify the proof of Theorem 5.1 in the following way. 

Let z = 1. Since {Eh)n,{uQ)h is strictly convex, for every ^ G d{Eh)r,Xuo)S'^h,i)i we have 

iEh)h,{uo)hiv) > {Eh)rj,(uo)hiul i) + {C,v- for every v G Eh- 

In particular, choosing v = {uo)h and recalling the definition of {Eh)r)^(uo)h we have 

Ehi{uo)h) > Eh{ul^^) + - {uo)h\£^ + iuo)h - VC G d{Eh)rj,(uo)A^h,i)- (7-5) 

Recall now that u\ i is the global minimizer of {Eh)n,(uo)h Therefore, there exists 

G Fh such that A^r^ G d{Eh)n,(uo)hi'‘^h,i)- Therefore, by ( [7^ , 

Ehi{uo)h) > Ehiui^i) + z?|4,i - {uo)h\£^ + {Alr\ {uo)h - 

At this point we can finally dehne the function qh^s ■ [0, T] —)• Fh- Since Uh^s is a discrete 
quasistatic evolution, for every z G N with id < T there exists g* G Fh such that G 
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d{Eh){uh, 5 {iS)). Then, we define qs : [0,r] —>■ Eh as 

te[0,6), 

f G [0, T] n [i5, {i + 1)(5), for i G N with i5 < T. 
Thus, repeating the proof of Theorem |5.1| for i = 1 we have 

- {uo)h\l^ + Eh{ul^i) - Eh{{uo)h) 

< {r^ - qs{0), - Ah{uo)h)Th + / A(s)>^h ds 

Jo 

= [ (/(O), A(s))j-^ds. 

Jo 

Recalling that Lemma |3.3| holds true also in this case, we have 

Eh{uh,5{S)) < Ehiui^i) < Eh{{uo)h) + / (q^s), fh{s))Th ds- 

Jo 


Since for i > 2 the proof of ( |5.l[ ) can be repeated with no modifications, this shows Step 1 of 
the proof of Theorem |5.1[ 


In order to prove condition (ii) of Theorem 5.1, we proceed in the following way. By the 
minimality property of u\i, we have 

Eh{ui^i) < Eh{ul^i) + r?|4,i - (uo)/i||^ < Eh{uJh{S)) + v\<^h{S) - {uo)h\£^ < C, 

where, by ( |7.1[ ) and ( 7.2| ), C is a constant depending only on rj, oj and no (and not on h and 6). 
Then, thanks to (7.11 

Eh{Uh,l) + \^hUh,l\Th ~ Eh{u\ i) + |/h(<^)l3^^ < C, 

where C is again a (possibly different) constant, depending only on rj, iv and uq (and not on h 
and 6). Thus, by condition (Jl) and by equicoercivity of the family of functionals [Eh)h^D 
have that 

< c 

for some constant independent of h. Therefore, thanks to property (J3) and recalling that 
Alr^ G d{Eh)r^,(^uo)^{ul^i) 

^ {Ehiul^i) + 1) + “ i^o)h\£n < C, 

for some constant C independent of h. From this, in particular, we obtain that 

k5(0)|j-^ = \r^\Th <C < -Eh{uh,s{0)) + c, 

which gives the equiboundedness of \qs{t)\j^h. ™ the interval [0,<5). At this point, the proof of 
Theorem |5.1| can be repeated without any modifications. 


Step 2: We apply Theorem 2.16[ 


By Step 1, properties (a)-(c) of Theorem 2.16 hold true. In particular, (a) implies (a’). By 

it 

(7.6) 


(b) of Theorem 2.16, there exists a bounded measurable function qh : (0, T] —>■ Eh such that 

A*qh{t) G dEh{uh{t)) 


for every t G (0, T\. 


Let now t G (0,T] be fixed. Thanks to Remark 2.2, we have 

for every w G ker(Aj^). 


0 < liminf 
£—>■ 0 + 


Eh{uh{t) + £w) - Eh{uh{t)) 
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A careful inspection of the proof of [6l Prop ositio n 3.1] shows that last inequality implies (b’). 
Let us now show (c’). From (c) of Theorem 


L^(0,T), and for every t £ [0,T] we have 


2.16 


the function s i—>■ {qh{s), fhis)) belongs to 


Ehiuhit)) < E{{uo)h) + / {qh{s),fhis))j^^ds. 

Jo 

Recalling that fh{s) = AhiOh{s) and that the linear operator Ah is independent of time, we have 

Eh{uh{t)) < E{{uo)h) + f {qh{s),fh{s))Th^s 

Jo 

rt 


= E{{uo)h)+ / {qh{s),AhU}h{s))Thds 
Jo 


= E{{uo)h) + / {Alqh{s),u}h{s))eh ds- 

Jo 


By (7.6), for every s £ (0,T) we have A*qh{s) £ dEh{uh{s)). Since Uh{s) £ for every 
s £ (0, T), by Remark 6.1 we have 


{J^hQhis),uJhis))e,, = / Vuhis) ■Vuhis)dx for every s G (0,r). 
Jn\r 


Therefore, 


Eh{uh{t)) < E{{uo)h) + / Vuh(s) • Va;h(s)dj;ds, 

Jo JQ\r 


which gives (c). Finally, property (d’) directly follows from Remark 2.17 and Remark 5.2 □ 

We can now pass to the limit as /i —)• 0''~. 


Proof of Theorem 7.5. Let t £ [0,r]. We will use an argument similar to the one in the proof 
of Theorem 12.161 

Step 1: Proof of (a) and (c). 

First of all, we fix the subsequence {hk)h^f^ given by 


hk ■- 


k£N, 


so that 


£h, C £h^ for every m> 1. 


(7.7) 


By (7.1), we have 


ujhf. ^ Jo strongly in L^([0, Tj; F7^(Q)) as A: —)• oo. 

Thus, there exists a set A2 C [0,r] with £^(A2) = 0 such that Lo{t) is well defined for every 
t £ [0, T] \ A2 and 

Johkit) uj{t) strongly in for every t £ [0,T] \ A2 as A: —)■ 00. (7.8) 

For every h £ D, let Uh ■ [0, T] —)• £h be given by Theorem |7.6[ We define 

[ Vuh^it) ■Vujh^^{t)dx for every t G [0,r] \ A2, 

Jn\r 

0 for every A G A 2 , 


6»fc(t) := < 
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and 


6{t) := limsup0fc(i) for every t G [0,T]. 

k^oo 


By definition of 0, for every t G [0, T] we can extract a subsequence (possibly depending 

on t) such that 

9{t) = lim 6k (t) for every t G [0,r]. 


J^OO 


By (7.4), for every t G [0,r] we can extract a further subsequence (not relabelled) such that 

Uh^. (t) u{t) weakly in \ F) as j —>■ oo. (7-9) 

for some u{t) G \ F) with ||u(t)||Hi(f^\r) — ^'2- By repeating what was done in the proof 

of Theorem 2.16, we can show that the subsequence can be chosen in such a way that 

the map u : [0,T] —)• ?7^(n \ F) is measurable, and this shows (a). 

Let us now show the energy inequality. By (7.8) and (7.9) we have that, for every f G [0, T] \A2, 

0(t) = limsup 0fc(f) = \im 9k {t) = lim / (f) • (f) dx = / Vu{t) ■Vu}{t)dx. 

fc^oo Jn\r ^ ^ Jn\r 

In order to prove that 9 G L^(0,T) we first observe that 9 is measurable, since it is the lim sup 
of a sequence of measurable functions. Moreover, we have 


rT fT 

\9{t)\dt= / 


n\F 
■T 


Vu{t) ■ Vuj{t) dx 


dt< ||Vu(t)||2,2(Q\r)l|Vw(t)|U2(Q\r) 


<C2 ||Va;(t)||i2(Q\r) dt < C2\/T’||w||L2((o,T);i4i(n))- 

Jo 


By (c’) of Theorem 7.6 we have, for every j G N and for every t G [0, T] 


(t)) = Eh. {uh. (t)) < E{uo,hk ^) + / / ^Uh. (s) ■ Vujh. (s) dxds. 

^ j •’Jo Jn\r ’ ’ 


(7.10) 


Note now that the energy E(-) is lower semicontinuos w.r.t. weak convergence in \ F). 

Moreover, since uo,hfe, strongly in \ F), we have 


lim E{uo,hk ) = E{uo). 


Therefore, taking the limsup in j of the (7.10), and using Fatou’s Lemma 


E{u{t)) < lim inf E{uh^ {t)) < F1(mo) + l™sup 


^--v-u; 1 —, , Vu/ij^,(s) • Va;/ij^,(s)dxds 

’ j^oo Jo Jn\r ^ ’ 



< E{uo) + limsup 



Vuh^{s) ■ '^iOh^is) dxds 

k—^oo Jo JQ\r 

< E{uo) + / limsup / Vuh^^is) ■VuJhk{s)dxds 
Jo k^oo JO\r 
ct 


— E(uq) + 


0 Jn\r 


Vu(s) ■ Va;(s) dx ds. 


so that (c) follows. 
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We finally prove the stationarity. Since Uh{0) = —)• mq as /i —)• 0, we only have to prove 

the condition at a point t G Let G H^{Q, \ L) with = 0 on Then (see for 

instance [Ml), we can find a sequence ) such that 


ieN 




tp strongly in H^{^} \ T) as j —>■ oo 


and iph,,. £ with = 0 on for every j G N. Note that, by (7.7), we have 
V'/ifc, £ ‘S’/ifc. with 'iphk, = 0 on do^ for every j > 1. 


Therefore, by (|7.3|) 

j 

ln\r 


[ Vixfe (t) • dx 

Jqw ^ ‘ 




dTi 


d-l 


for every j > 1. By (7.9) we have 


lim / dx= / Vu{t) ■V'lph^.dx. 

Jn\r ^ ' Jn\r ' 

Define now, for every t G [0,r] and for every j > I, the function fj{t) : T —; 
fj{t) := -[V'h,j5'(|[wftfc.(i)]|)sign(K,,(i)])lj„^ W “ 


(7.11) 


(7.12) 


as 




We want to prove that for every t G [0, T] 

Up to extracting a further subsequence, we can assume that 

liminf/j(t) = lim fj{t) Ti^-a.e. in T, 
j^oo j—^oo 


Ti^-a.e. in T. 


and 


lim [uh. (t)] = [u(t)] 


"H^-a.e. in T. 


j->oo 


(7.13) 

(7.14) 

(7.15) 


Now, let us fix X G Ju{t) such that (7.14) and (7.15) hold true. Then, for j G N large enough we 
have 

X e Ju. (t) and sign{[uh^Xt)]{x)) = sign([u(t)](x)). 

fcj J 

Therefore, 

liminf/j(t)(x) = lim fj{t){x) 


J^OO J—^OO 




= -[V’fefcJ(4:)5''(|[w(i)](4;)|)sign([u(t)](x))lj„(^^(x) - |[V’ft^J(x)|ljc^^j(x 


(7.16) 


for H^-a.e. x G T n Ju{t)- If, instead, x G then recalling that 0 < < 1 we have 

liminf/,(t)(x) = lim fj{t){x) 

j^OO 

= lim ](x)(/'(|K ,(t)](x)|)sign([uft .(^)](®))lu^ (0(2^) - l[V’/ifcJ(a:^)|lj= (x) 

j^oo ‘ J J i^kA ' ' ‘ hi. , v; 


> -|[V’/^feJ(a:)| = -liXhk^KxMlj- (x). 


(7.17) 
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Combini ng (|7.16[ ) and ( |7.17[ ) we obtain ( |7.13[ ). Thanks to ( |7.12[ ) and ( |7.13[ ) we can pass to the 
limit in (|7.11|), obtaining 


f Vu{t) ■ V'lphk, dx = lim f Vuhf^Xt) ■ V'lph. dx 
Jn\r ' Jn\r ^ ' 

> [ liminf f-[V’ftfe,]ff'(|K,^(t)]|)sign([n/,,Xi)])lj„,, m “ (,) ) ^77 


/p j->-oo 


Finally, passing to the limit as Z —>• oo we have 


d-l 


Jn\r 
and we conclude. 


Vu(t) • VV’dx > (|[^^(i)]|)sign([u(t)])lj„(,^ 


, 1 dTi 

u(t) 


d-l 


□ 


8. Numerical experiments 


The scope of this section is to practically show that the procedure illustrated in the previous 
sections can be effectively implemented and produces the desired quasistatic evolution, according 
to the one described in [6j. We refer the reader to Section]^ for the notations used here. 

8.1. Numerical simulations in 1 dimension. We first analyze the results obtained for a one¬ 
dimensional problem, when fl C M. Despite its simplicity, the one-dimensional setting allows us 
to give a detailed comparison between numerical results and analytic predictions, since in this 


case the explicit solutions of (6.5) are known. We consider the following geometry: 

n=[0,27], 7 = 0.5, r = {7}, dDn = {0,2£}. 

We follow the evolution in the time interval [0, T] = [0,1], and the external load applied to the 
endpoints = {0, 27} is given by 

uj{t){x) = 2{x — i)t, for every x G [0, 27] and t G [0,1]. 

We uniformly discretize the domain into 2N = 80 intervals, so that the spatial discretization 
step is given hy h = i/N. Finally, we choose a time step S = 0.02, so that the total evolution is 
concluded after 50 time steps. In our specific case, by Proposition |6.3| and a direct calculation 
we have that the parameter rj in condition (J2) can be taken as 

V = ^ +max{4,4\/7}. 


where the constant R is the one appearing in the definition of the function g, see (6.2). 


From a practical viewpoint, the computational time needed to solve the minimization problem 


(3.3) could grow without any control. Hence, for any i G {0,... ,50} and j G Nq fixed, we stop 
the minimization loop as soon as ||Hu —/(i(5)|| < 10~®, where A is the trace operator and / = Aui 
(we omit the dependence of H on /i to ease notation). That is, Vj is chosen in such a way that 
||Hu* — f{iS)\\ < 10“®. Instead, we stop the external loop (that is, the limit of u* as j —>• oo), as 


soon as 


\VJ-V 


i-il 


<10 The main reason for these choices is that the quasistatic evolution 
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generated by the algorithm is extremely sensitive to any perturbation. Thus, a larger time step 
5, or a too badly approximated critical point at each time step, could lead to nonphysical results. 

We observe that the evolutions discussed in [6l Section 9] depend on the size of the parameter 
R. Therefore, in order to compare our results with those in [6], we distinguish two cases. 

Case R > 21. When R is chosen large with respect to the size 21 of the elastic body, the 
evolution found by numerical simulations evolves along global minimizers of the energy, and we 
can observe the three phases of the cohesive fracture formation: non-fractured, pre-fractured 
(that is when the opening of the crack is smaller than R and cohesive forces appear), and 
completely fractured (when the opening of the crack is larger than R and the cohesive forces 
disappear), see Figure]^ Note that in the time interval [0,0.5] the evolution follows the elastic 



Displacement at t = 1 

1 . 51 -^^^- 


0.5 

0 

- 0.5 


0 0.2 0.4 0.6 0.8 1 


Figure 4. The evolution of the quasistatic cohesive fracture for R> 2i at time 
instances t = 0, 0.5,0.52,1. 

deformation. After t = 0.5 a fracture appears, since the elastic deformation is not any more a 
critical point of the energy functional (see [H Section 9]). Then, the pre-fracture phase starts, 
showing a bridging force acting on the two lips of the crack. At time t = 1 the cohesive energy 
reaches its maximum, and the body is completely fractured. It is worth observing that this 
evolution coincides with the one analytically calculated in [6l Section 9]. We can also investigate 
what happens from the energy point of view, see Figure We have a smooth transition between 
the different phases, and the total energy has a nondecresing profile. The beginning of the pre¬ 
fractured phase can be observed at the 25*^ time step (i.e. at time t = 0.5), when the elastic 
energy (in red) starts decreasing and the crack energy (in blue) starts increasing. The final 
phase of complete rupture is then attained at the final time step t = 1. Although we focused on 
the time interval [0,1], one could check that the three energy profiles remain constant for t > 1. 

Case R < 2i. The evolution of the system changes radically when R < 2L In this case 
(see Figure]^ the failure happens instantaneously, without a bridging phase, and thus the body 
exhibits what in literature is known as brittle behavior. More precisely, in the time interval 
[0,0.5] the evolution follows again the elastic deformation, and a crack appears at t = 0.5. 
However, immediately after t = 0.5 the body is completely fractured, and no cohesive forces 
appear. It is important to observe that in this case we actually observe an evolution along 
critical points that are not global minimizers. Indeed, the evolution is elastic until t = 0.5, 
although it would be energetically convenient to completely break the body at some earlier time 
t < 0.5 (see [H Section 9] for a detailed description of all critical points). Thus, we see that 
the algorithm chooses the critical point which is the closest to the initial configuration, even if 
other options are available, which are more convenient from an energetic point of view. This 
evolution is particularly supported by the idea that in nature a body does not completely change 
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Figure 5. The total, fracture, and elastic energy evolution of the quasistatic 
cohesive fracture for R > 2i. 



Figure 6. The evolution of the quasistatic cohesive fracture for R < 2£ at time 
instances t = 0, 0.5,0.52,1. 

its configuration crossing high energetic barriers if a stable configuration can be found with less 
energetic effort. 

Also in this case, we can observe the evolution from the energetic viewpoint, see Figure [7| At 
time t = 0.5, when the elastic deformation ceases to be a critical point, the domain breaks and 
the total energy decreases up to the value of R/2, so that no bridging force is keeping the two 
lips together. As we already observed, the evolution along global minimizers would instead lead 
to a fracture way before the critical load is reached. 

Again, the evolution found with our numerical simulation coincides with that one given in 
[U Section 9]. In particular, our simulations agree with the crack initiation criterion (see [6l 
Theorem 4.6]), which states that a crack appears only when the maximum sustainable stress 
along F is reached. In this case, this happens at t = 0.5, when the slope of the elastic evolution 
reaches the value = 1. 

8.2. Numerical simulations in 2 dimensions. Having a first analytical validation of the 
numerical minimization procedure, we can now challenge the algorithm in the simulation of two 
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Figure 7. The total, fracture, and elastic energy evolution of the quasistatic 
cohesive fracture for R < 2i. 


dimensional evolutions. We now consider the domain introduced in Section 6.2 setting I = 0.5, 
2N = 8, and n = 1/2. Within this choice, the crack initiation time is reduced exactly of a 
factor 1/2, allowing us to speed up the failure process. Since all the computations are performed 
on a MacBook Pro equipped with a 2.6GHz Intel Core i7 processor, 8GB of RAM, 1600MHz 
DDRS, the two dimensional simulations are performed only for a qualitative purpose. Indeed, 
we are mainly interested in showing that our algorithm produces physically sound evolutions 
also in dimension 2, and when the external displacement / is non-trivial. The very sparse 
discretization of the domain D is due to the fact that the minimization in (3.3) requires a huge 
computational effort, both in terms of time and memory. Indeed, in order to implement more 
realistic experiments, with a finer discretization, we would need to modify the architecture of 
the minimization algorithm, in such a way that it may run on parallel cores. 

We perform two different series of experiments, one with boundary datum 


t(Ji(t)(x) = 2(xi — tjt, for every t G [0,1] and x G D, 
see Figure and the other one with boundary datum 

uj 2 {t){x.) = 2tcos ^ (xi — £), for every t G [0,1] and x G D, 

see Figure]^ Here, we denoted by x = (xi,X 2 ) the generic point of D = (0,1) x (0,1). We now 
need to reduce the tolerance of the termination condition of the outer loop of the Algorithm, 
setting it to 5 • 10~^^. Indeed, we experimented that for bigger values of this tolerance some 
instabilities in the solution were introduced, leading to an asymmetric evolution, also in the case 
of oji as external displacement, where we expect an invariant behavior with respect to the space 
variable X 2 - 


Case R > 21. In Figure and we report 4 different instances of the evolution for the two 
different boundary data, when R > 2i. When the external displacement is wi, which is constant 
with respect to the second coordinate X2, we observe that the evolution is also constant with 
respect to X 2 - For both boundary data, the failure of the body undergoes the three phases of 
deformation, as it happened in the one dimensional case. 
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Figure 8. The evolution of the quasistatic cohesive fracture for R>2^aX time 
instances t = 0,0.24,0.32,0.5 with external displacement wi. 



Figure 9. The evolution of the quasistatic cohesive fracture for R>2^sX, time 
instances t = 0.1,0.24,0.34,0.5 with external displacement U 02 - 


Case R < 2i. When the boundary datum is wi, see Figure [Tol the specimen breaks in a 
brittle fashion, without showing any cohesive intermediate phase. This simulation is actually an 
evidence that the algorithm still characterizes the correct critical points, following the principle 
that the domain should not fracture as long as a non-fractured configuration is still a critical 
point. We conclude commenting the simulation where the boundary datum is 012 with R < 2£, 
see Figure [m By setting a displacement highly varying with respect to the X 2 coordinate, we 
observe that the different phases of the fracture formation can cohexist. At time t = 0.24 the 
domain still presents no fracture, as expected by the previous numerical experiments. Then, at 
t = 0.34, a pre-fracture appeares, but only at those points where the external load is bigger, i.e. 
around X 2 = In fact, even at the final time t = 1, the domain is not completely fractured. 

Note that, when the boundary datum is a;i, the evolution coincides with the one obtained 
analytically [U Section 9]. In particular, the fracture appears at t = 0.25, when the slope of 
the elastic evolution reaches the value Kg'{0) = 1/2 and thus the crack initiation criterion is 
satisfied. 




Figure 10. The evolution of the quasistatic cohesive fracture for R < 2t at time 
instances t = 0.04,0.24,0.26,0.5 with external displacement wi. 
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Figure 11. The evolution of the quasistatic cohesive fracture for R <2^ time 
instances t = 0.02,0.22,0.24,0.5 with external displacement u) 2 - 
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